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Q_i' Abstract: 

(D ■ 

'^ ' We develop a new theoretical framework for the description of leading twist light-cone 

baryon distribution amplitudes which is based on integrability of the helicity A = 3/2 
/\ ■ evolution equation to leading logarithmic accuracy. A physical interpretation is that 

d • one can identify a new 'hidden' quantum number which distinguishes components in the 

A = 3/2 distribution amplitudes with different scale dependence. The solution of the 
corresponding evolution equation is reduced to a simple three-term recurrence relation. 
The exact analytic solution is found for the component with the lowest anomalous di- 
mension for all moments N, and the WKB-type expansion is constructed for other levels, 
which becomes asymptotically exact at large A^. Evolution equations for the A = 1/2 
distribution amplitudes (e.g. for the nucleon) are studied as well. We find that the 
two lowest anomalous dimensions for the A = 1/2 operators (one for each parity) are 
separated from the rest of the spectrum by a finite 'mass gap'. These special states can 
be interpreted as scalar diquarks. 
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1 Introduction 

There exists a general consensus that exclusive processes involving large momentum 
transfers are dominated by 'valence' components in hadron wave functions with the 
minimum number of Fock constituents [0, 3. It is equally generally accepted that the 
asymptotic behavior of exclusive amplitudes is in most cases determined by the so- 
called 'hard-rescattering' mechanism involving configurations of partons with almost 
zero transverse separations, although the theoretical status of the dominance of small 
transverse distances is somewhat weaker for baryons [^ than for pions [^ Q. 

As always in a field theory, extraction of the asymptotic behavior introduces diver- 
gences. In the present context, infrared divergences in perturbative diagrams describing 
the hard rescattering are removed by renormalization of nonperturbative scale-dependent 
distribution amplitudes which are defined in terms of the Bethe-Salpeter wave functions 
integrating out transverse degrees of freedom 

\ki,±\<tJ- 

0(xi,...,x„;/i) ~ / d'^ki^±...d'^kn,±^Bs{xi,kix,---;Xn,kn,±) (1.1) 

with Xi being the longitudinal momentum fractions carried by partons. The concept 
of distribution amplitudes is central for the theory of hard exclusive processes where 
their role is analogous to that of more familiar parton distributions in the description of 
inclusive processes. 

The theoretical basis for studies of distribution amplitudes [^] is provided by their 
definition in terms of hadron-to- vacuum transition matrix elements of non-local gauge- 
invariant light-cone operators of the type (in a suitably chosen gauge) 

q{zin)q{z2n), 

eijkq\zin)q\z2n)q^{zzn) (1.2) 

for meson and baryon distributions, respectively. Here g* is a generic quark field with 
the color i, n is an auxiliary light-like vector n^ = and Zi are real numbers that 
specify quark (antiquark) separations. More specific definitions will be given below. 
The nonlocal operators as above are understood as generating functionals for the series 
of local operators obtained by their Tailor expansion at short distances (contraction of 
the derivatives with the light-like vector ensures taking the leading twist part) and the 
precise relation is such that moments of distribution amplitudes are given by matrix 
elements of the contributing local operators [Q. The scale dependence of the moments 
of distribution amplitudes corresponds to the renormalization group (RG) evolution of 
local operators and can be studied using familiar methods. 

The specific problem for distribution amplitudes is to take into account the addi- 
tional mixing with operators containing total derivatives that cannot be neglected in 
contrast to inclusive processes where only forward matrix elements are being consid- 
ered. It was noticed that the RG evolution is driven to leading logarithmic accuracy 



by tree- level counterterins which thus have the symmetry of the bare QCD Lagrangian 
and, in particular, the conformal symmetry. As a consequence, operators belonging to 
different irreducible representations of the conformal group cannot mix under renormal- 
ization in one loop @, |^, 0, P, H- This observation solves the mixing problem for meson 
(two-particle) distributions in which case a single independent local conformal opera- 
tor exists for each moment. The corresponding anomalous dimension can be continued 
analytically to non-integer (complex) moments, defining the Altarelli-Parisi evolution 
kernel: coefficients in the expansion of meson distributions in the basis of Gegenbauer 
polynomials are renormalized multiplicatively and with the same anomalous dimensions 
as in deep inelastic scattering [Q. Consequently, assuming 'good' behavior at complex 
infinities, the distribution amplitude can be restored by inverse Mellin transform from 
analytically continued values of the moments; hence the partonic interpretation of the 
distribution amplitude proves to be consistent with its renormalization properties (scale 
dependence) . 

The three-quark baryon distribution amplitudes bring in a complication of principle. 
The conformal symmetry allows to resolve the mixing with total derivatives but is not 
sufficient to diagonalize the mixing matrix completely. For fixed operator dimension, 
alias for fixed total number A^ of covariant derivatives, D^ = d^ — igA^, there exist A^-|- 1 
independent local operators (modulo operators with the total derivatives) 

ON,k = q (d -nfq {n- of-'q, k = 0,...,N (1.3) 

corresponding to A^ + 1 genuine independent degrees of freedom of the three-quark sys- 
tem. One is left with a nontrivial (A^ + 1) x (A^ + 1) mixing matrix that has to be 
diagonalized explicitly order by order; see, e.g., [|iy, 0, [T^, |, |r^, |I^. The resulting 
A^ + 1 multiplicatively renormalizable operators have different (in general) anomalous 
dimensions whose analytic expressions are not known. Apart from mathematical incom- 
pleteness, absence of analytic results means that the general structure of the spectrum 
is unknown and, in particular, analytic continuation of the anomalous dimensions to 
complex moments A^ is not possible. This, in turn, implies that partonic interpretation 
of different 'components' in baryons is not understood beyond the tree level. 

This problem was well known but considered as a relatively minor one and did not at- 
tract due attention in the past. One reason was that the scale dependence of distribution 
amplitudes turned out to be rather mild in a perturbative domain and it seemed prema- 
ture to elaborate on the evolution before gross features of nonperturbative distributions 
were understood at low scales. We think that this logic is flawed and the general experi- 
ence with hard processes in QCD rather suggests that 'intrinsic' parton distributions at 
scales of order 1 GeV cannot be viewed as purely nonperturbative and disconnected from 
perturbative evolution. Despite an obvious fact that perturbative calculations cannot 
be made quantitative at low scales, there is increasing evidence that general patterns of 
the perturbative gluon emission are continued to very low momenta. For example, the 
shape of deep inelastic structure functions at 1 GeV appears to be largely determined 
by perturbative soft gluon radiation. Small differences in the perturbative evolution of 



different components in nucleon distribution amplitudes are strongly amplified in the 
nonperturbative domain and one may think that such differences build up gross features 
of distribution amplitudes at scales of order 1 GeV, where from the perturbation theory 
becomes quantitative. Viewed from this perspective, a detailed study of the evolution of 
baryon distribution amplitudes becomes mandatory. 

In this paper we suggest a new approach to the construction of baryon distribution 
amplitudes which is based on the recent finding ||T5| that the evolution equation for the 
baryon distribution amplitudes with maximum helicity A = 3/2 is completely integrable. 
That is it possesses a nontrivial integral of motion which we identify as a new 'hidden' 
quantum number that distinguishes components in the A = 3/2 distribution amplitudes 
with different scale dependence. It is interesting to note that the A = 3/2 evolution equa- 
tion is equivalent to the quantum mechanical problem that has already been encountered 
in QCD in the studies of the Regge asymptotics of the scattering amplitudes [jl6|, |17 



and in the theory of integrable models as the so-called Heisenberg XXXs=_i spin magnet 



||T6| , |T8|. This problem has been studied in some detail using nontrivial mathematical 
methods and the results can be adapted to the present context. 

Our approach is advantageous compared to the standard formulation in at least 
two aspects. First, from practical point of view an important simplification is that 
diagonalization of a (A^ + 1) x (A^ + 1) matrix is replaced by solution of a simple three- 
term recurrence relation, which reduces the computer time significantly. Second, we 
obtain explicit analytic solutions to the evolution equations in all important limits. In 
particular, we will be able to identify trajectories of the anomalous dimensions and 
calculate them (and the corresponding eigenfunctions) using a WKB type expansion for 
large values of A^. 

These results apply in full to the A = 3/2 distribution function of the A-resonance 
and allow for a fairly complete description. The evolution equation for the A = 1/2 
distributions is not exactly solvable, but the difference to the A = 3/2 evolution can 
be considered as a small (calculable) perturbation for the most part of the spectrum. 
On the other hand, the structure of the lowest eigenstates is changed drastically. As 
we will demonstrate, the two lowest anomalous dimensions for the A = 1/2 operators 
decouple from the rest of the spectrum and are separated from it by a finite 'mass gap'. 
These two special states (one for each parity) can be interpreted as bound states in the 
corresponding quantum-mechanical model, and, somewhat imprecisely, can be thought 
of as corresponding to formation of scalar diquarks. 

As a byproduct of our study, we construct a convenient orthonormal basis for the 
expansion of three-particle distribution amplitudes, which is more suitable compared to 
standard Appell polynomials. 

The presentation is organized as follows. Section 2 introduces definitions and the 
general framework for the construction of baryon distribution functions and their renor- 
malization. Section 3 is devoted to the conformal symmetry of the evolution equation 
and the conformal expansion of distribution amplitudes. Section 4 presents a detailed 
study of the exactly solvable evolution equation for the maximum helicity A = 3/2. In 



Section 5 we consider the evolution equation for the A = 1/2 distributions. A short 
summary of main results of phenomenological relevance is given in Section 6 and the 
general conclusions in Section 7. Appendix A presents an explicit construction for the 
Racah 6j— coefficients of the SL{2) group and in Appendix B we consider the effective 
Hamiltonian for low-frequency modes of the A = 1/2 evolution equation. 



2 General framework 
2.1 Distribution amplitudes 



Following Refs. |r^ we define the leading twist nucleon distribution amplitude as the 
corresponding matrix element of the gauge-invariant three-quark nonlocal operator 



{0\ui{zi)u^^{z2)d'^{z3) Ui'iizi, zo)Ufj{z2, zo)Uk'k{z3, zo)e'^^\P{p, A)) 



(2.1) 



X {i^C)apil6N),Viz,p) + {^^,C)^(,N,Aiz,p) + (ta^^p^Q^p {^^^,N),T{z,p) 
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where a^^u = f [7/^,71/], C is the charge conjugation matrix, \P{p,X)) is the proton state 
with momentum p and helicity A, and A^ is the proton spinor. All the interquark separa- 
tions are assumed to be light-like, e.g. u{zi) denotes the u-quark field at the space point 
Zin with n^ = 0, and U{zn, Zq) are non-Abelian phase factors (light-like Wilson lines) 



U{zn, Zq) = P exp 



ig / dt {zn - zo) n^A^'{tZn + (1 - t)^o) 



(2.2) 



Because of the light-cone kinematics, the matrix element in fact does not depend on zq 
and the phase factors can be eliminated by choosing a suitable gauge. To save space we 
do not show the gauge phase factors in what follows, but imply that they are always 
present. 

The invariant functions V, A, T have the following symmetry properties []19| 



r(l, 2, 3) = V{2, 1, 3), Ail, 2, 3) = -A{2, 1, 3), r(l, 2, 3) = r(2, 1, 3). (2.3) 
and can be expressed in terms of a single function 07v as 

2T(1, 2, 3) = 0^(1, 3, 2) + 0^(2, 3,1), 

0^(1,2,3) = V(1,2,3)-A(1,2,3). (2.4) 

Here (J)n is the leading twist proton distribution amplitude. If it is presented in the form 

^„(..rt./p.exp[-,E...(p-«)]^«(..). (2.5) 

where 

(2.6) 



Vx = / dxidx2dx3 6{l — Xi — X2 — x^) 



then the variables Xi have the meaning of the longitudinal momentum fractions carried 
by the three quarks in the nucleon, < Xj < 1 and ^Xj = 1. 

For what follows, it is convenient to introduce quark fields with definite chirality 

gW) = i(l±^5)g. (2.7) 



The definition in (|2.1| ) is equivalent to the following form of the proton state 0, ^ 



\P{p, A = +1/2)) = hj ^^^/^"^^^^^^ { \uHx,)u^{x,)d\xs)) - \u\x^)d^ix2Wixs))}, 

(2.8) 
where the standard relativistic normalization for the states and Dirac spinors is im- 
plied [pT| . The distribution amplitude (pN can be defined in terms of chiral fields: 



= -^fN{pn) fiN\p) Vx exp[- ipn{zixi + Z2X2 + Z3X3)](j)N{xi, fi"^) (2.9) 
so that moments of (pN 

(f>N{h,k2,h) = / 'E'xXi^X2^X3^ 0Ar(xi,X2,X3,/i^) (2.10) 



can be calculated as reduced matrix elements of the renormalized three-quark leading 
twist operators 

Ol[\,M = {nDf'u\Q){Cfi){nDf-u^mnDf^ 0{S)), (2.11) 

where D^ = d^ — igA^ is a covariant derivative. 

The leading twist distribution amplitudes of the A resonance can be obtained in the 



similar manner ||2^. For definiteness, we will consider the distribution amplitudes of 
A"'"^ only, all the other ones can be reconstructed with the help of the isospin symmetry, 
see [EDI. For this case one writes 



(0|<(^0«^,(^2)tx?(^3)|A++(p,A))6^^-^ = (2.12) 

,1/2 

- \ {io,,C)^^ (7,A'^),Ta(z,p)| - )jf iia^^C)^^ [p,^^ - \m^^,N^ <^T{z,p). 

Here A^(j9) is the A resonance spin— | vector: 

(^-Ma)A^ = 0, A^A^ = -2Ma, 7/.A'^(p)=PmA'^(p) = 0. (2.13) 



The dimensionless amplitudes Va{x),Aa{x),T/^{x) determine the distribution of quarks 
in the A(|A| = 1/2) state and satisfy the following symmetry relations 

Va(1, 2, 3) = Va(2, 1,3), Aa(1, 2, 3) = -Aa(2, 1, 3), 

Ta(1,2,3)=Ta(2,1,3), (2.14) 

TA(l,2,3) = yA(2,3,l)-v4A(2,3,l). 

Therefore, only one function is independent; we choose (cf. pO| ) 

1 /2 
0A {X1,X2,X3) = Va{xi,X2,X3) - Aa{xi,X2,X3) =Ta{x3,Xi,X2) (2.15) 

as the distribution amplitude of the A(|A| = 1/2) resonance. The remaining function 

3/2 

0^ is totally symmetric in all its arguments and determines the distribution of quarks 
in the A(|A| = 3/2) state. 

The structure, again, becomes more transparent when going over to chiral quark 
fields. The definition in Eq. ( p.l3| ) is equivalent to the following structure of the A 



resonance states pO| 



a/2, 



\AiX = l/2)) = fj' Ivx^^^t^ju\x,y{x2)u\x3)), 

|A(A = 3/2)) = fTl'^^-^^^==\^'(^^W(^2)u\x3)), (2.16) 



p3/2 _ /2 xl/2 



where /^ = a/|A^ /Ma, and the distribution amplitudes can be defined through the 
nonlocal matrix elements: 

= --XAn''{fiA)l / Vx exp [ - ipn{ziXi + Z2X2 + Z3X3)] (pAixi, /i^) (2.17) 

and 

{0\e^^'{u]{z,)Ca,,n,u]{z2)){A''ji)ul{z3)\A{p)) = 

= i /a (P"-)^! ^A^''' / Vx exp [ - ipn{zixi + ^2X2 + Z3X3)] 0a (2^*) A^^)- (2-18) 

2.2 Renormalization 

In this paper we will be interested in the scale dependence of the baryon distribution 
amplitudes. For each of them, = 0^, 0^ , 0^ , we anticipate an expansion of the type 

(t){Xi,lI^)=XiX2X3y^(j)nPn{Xi)\^ T , (2.19) 

^ V"s(/io)/ 



where 60 = ll/3A^c — 2/3n/, -P„(xj) are certain polynomials, 7„ are the corresponding 
anomalous dimensions and (pnifJ'o) are dimensionless nonperturbative parameters. The 
prefactor xia;2a;3 suggests the vanishing of the distribution amplitude at the end points 
Xk = and as we will show its presence is closely related to the conformal invariance 
of the evolution equations. Finding 7„ and P„ corresponds to explicit diagonalization of 
the mixing matrix of the three-quark composite operators and is fully equivalent to the 
solution of the corresponding Brodsky-Lepage equations 0. 

Renormalization properties of the relevant three-quark operators are most conve- 
niently presented in terms of their generating functionals (nonlocal operators) with three 
spinor indices |p2|| : 

B%\iz,, Z2, zs) = e'^\JiqjUz^n)i^q})^iz2n)iJiqlUzsn), (2.20) 

Bl%iz„ z,, Zs) = e''\jiqJUz,n)ijiq^^)piz2n)ijiqlUzsn), (2.21) 

with qi being a quark field of color i. Bs/2 gives rise to the distribution amplitude 
0^ and Bi/2 is relevant both for 0^ and 0^ . The nonlocal operators -B3/2 and 
Bi/2 do not mix with each other since they belong to different representations of the 
Lorentz group: (3/2,0) and (1,1/2), respectively^. For most of the discussion we will 
assume that all three quarks q have different flavor. Identity of the quarks does not 
influence renormalization but rather introduces certain selection rules which pick up the 
eigenstates with particular symmetry, to be detailed later. 

The renormalization group equation for the nonlocal operators ( |2.20| ), ( |2.21| ) can be 



written as 23, 22 



where EI is some integral operator corresponding, to the one-loop accuracy, to contribu- 
tions of the Feynman diagrams shown in Fig. |l]. 

To simplify notations, we factor out the color factors and trivial contributions of the 
self-energy insertions: 

e = (1 + 1/Nc)n + 3Cf/2, (2.23) 

with Cp = (A^c " l)/(2A^c)- It is easy to see that the gluon exchange diagram in Fig. [po 
vanishes unless the participating quarks have opposite chirality. The renormalization 
of the A = 3/2 operator Bs/2 ( |2.20D is therefore determined by the vertex correction in 



Fig. |l]a alone (in Feynman gauge). By explicit calculation one finds [|1^, ^, [Tl|, |T^, |T3[ 

^^3/2 = ^12 + ^23 + ^13) (2.24) 

where TiJ'^ are the two-particle kernels involving the i-th and k-th quarks, for example, 

ni,Biz,) = - [ —\a[B{z'^„Z2,Zs) - B{z,,Z2,zs)] 
Jo a I 



^ The transformation properties can be made manifest by going over to the two-component spinors 
and 7-matrices in the Weyl representation U^. 
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imfytyiT 



b 



Figure 1: Examples of a 'vertex' correction (a), 'excliange' diagram (b) and self-energy in- 
sertion (c) contributing to the renormalization of three-quark operators in Feynman gauge. 
Path-ordered gauge factors are shown by the dashed lines. The set of all diagrams includes 
possible permutations. 



+ a[B{zi,Z2i,Z3) -B{zi,Z2,Z3)]j, 



(2.25) 



with a = 1 — a and z'^ 



ik 



ZiU + z^a. 



In the case of S1/2 the vertex correction remains the same, but one has to add 
contributions of gluon exchange between the quarks with opposite chirahty. One obtains 



n 



1/2 



/L3/2 — rii2 — rL 



23' 



(2.26) 



where we assume that the first and the third quark have the same chirahty, as in (|2.21 



The kernels Hfj^ act on i-th and A;-th arguments of the nonlocal operators only, and can 
be written in the form 



ni2B{z,) 



VaB{z'^^,z^l,Z3) 



(2.27) 



with the integration measure Va defined in ( |2.6| ). 

Going over to local operators corresponds to the Taylor expansion of the generating 
functionals at small distances: 






B{z,,Z2,z,) = J2 E -r^f^VJ^nD)'''imnDr^qmnD)'^-^q{0). (2.28) 



N ki+k2+k3=N 



hi /C2! /C3 



The total number of derivatives N is preserved by the evolution so that the integro- 
differential equation ( |2.22| ) takes the matrix form, with the square matrix of size A^(A^ + 
l)/2 for each given N subsector. 

A generic local operator with A^ derivatives can be written as sum of monomials 
entering the expansion (|2.28| ) with arbitrary coefficients 



0= Yl c,,M2M{nDf^m{nDf'qmnDf^q{Q), 

ki+k2+k3=N 



(2.29) 



and can be represented by a polynomial in three variables 

^(xi,a;2,X3) = Yl Ck,k2kA^A^xl\ (2.30) 

fcl+A;2+A:3=Af 

In what follows we refer to '^{x.i) as coefficient function of a local operator. To justify 
the name, note that \E'(xj) serves as a projector separating the contribution of the local 
operator Oi^, to the nonlocal operator B{zi), which can be made explicit by writing 

O^ = ^idi,d2,d3)B{zuZ2,Z3)U^=o. (2.31) 

Local operators having the same number of derivatives all mix together so that the 
size of the mixing matrix for given N is A^(A^ + 1)/2. Since a local operator is completely 
determined by its coefficient function, diagonalization of the mixing matrix for operators 
can be reformulated as diagonalization of the mixing matrix for the coefficient functions. 
Requiring that (9^ ( |2.31|) is multiplicatively renormalized, one ends up with a matrix 



equation in the space of homogeneous polynomials of degree N of three variables^ 

n ■ ^N,q = SN,g'^N,g , (2.32) 

whose eigenvalues correspond to the anomalous dimensions 

7iv,g = (1 + 1/iVc) £N,g + 3/2 Cp. (2.33) 

Note that the eigenfunctions and the eigenvalues have two indices: A^ which refers to 
the degree of polynomial alias the total number of derivatives, and q which enumerates 
the energy levels. In the case of 7^3/2 we will later identify q with a conserved charge. 
The SL{2) symmetry of the equation ( p.32|) (see below) implies that the anomalous 
dimensions take real quantized values and the corresponding eigenfunctions are mutually 
orthogonal with the weight function a;iX2X3 

/ r'xXiX2X3 '^N,q{Xi)^!N,q'{Xi) ~ 6q^qi . (2.34) 

The same property allows to identify the eigenfunctions as the polynomials entering the 
expansion of the distribution amplitudes in ( p.l9| ) 

Pn{Xi) = ^NA^^). (2.35) 

To see this, consider the matrix element (0|C*|i?) = \l'(9i, (92, (93)(0|i?(zi, Z2, ^3)|-B)|z,=o 
and use the representation similar to ( P^.18D to get 

{0\O^\Bip)) ~ Jvx [v]/(9i,92,53)e-P"(^i-i+^^-^+^^-3)]^^^^ 0b(x„/x2) 

= {-ipn)^ Vx'^{xi,X2,X3)(f)B{xi,f/). (2.36) 



^Notice that the action of the evohition kernel H on the space of the coefficient functions is different 
from that on the nonlocal operator B{zi) 
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Substituting the expansion ( |2.19| ) into this relation and taking into account that the 
operator 0,i, is renormahzed muhiphcatively (by construction), one immediately finds 
that, first, the polynomials P„ have to coincide with '^N,q up to arbitrary normalization 
and, second, the nonperturbative coefficient (pn is given by the reduced matrix element 
of the operator 0^. 

Note that the 'Hamiltonians' 7^3/2 and 7ii/2 acting in the space of coefficient func- 
tions in ( |2.32| ) are not the same as those acting on nonlocal operators, although they are 
related, of course, and the precise connection can easily be established. Explicit expres- 
sions for 7^3/2 and 7ii/2 in the matrix representation can be found in |jlO|, |12|, ||, |13|. 

3 Conformal invariance 

The Lagrangian of massless QCD is known to be invariant under conformal transforma- 
tions. This symmetry survives for evolution equations at one-loop level since breaking of 
the conformal Ward identities induced by the nonzero trace of the stress-energy tensor 
is proportional to the QCD /5 function and is of order a^ B HI- One should expect, 
therefore, that the evolution operator 7i introduced in the previous section has the same 
symmetry and in particular commutes with the generators of the conformal group. This 
property imposes strong constraints on a possible form of the eigenfunctions: In a generic 
situation the eigenfunctions of two-particle operators are uniquely determined by confor- 
mal invariance whereas for three-particle operators one is left with an arbitrary function 
of one variable. Aim of this section is to work out the necessary framework. 

3.1 Collinear subgroup SL{2,M.) of the conformal group 

Algebra of the full conformal group contains the generators of dilatations D and special 
conformal transformations K^ in addition to the Poincare generators P^ and M^i,. The 
algebra reads 

[D, K^] = tK^, [K^, PJ = -2t{gf,,B + M^,), [D, M^,] = 0, 

[D, P^] = -zP^, [K^, M,,J = tig^^K, - g^^K^), [K^, K,] = 0, (3.1) 

plus usual relations for Poincare generators. Action of these generators on an arbitrary 
quantum field $ (e.g. quark or gluon) is given by (see, e.g. pi, M) 



[P^, <l>(^)] = -td^^z), [M^, , <l>(^)] = [t{z,d^ - z^d,) - S,,,] <l>(z), 

[D,<l>(z)] = -z(^^9, + /)$(^), 

[K^, ^{z)\ = -I {2z^z''d, - z^d^ + 2lz^ - 2zz^S^,) $(2) . (3.2) 

Here / is the canonical dimension of $ (/ = 3/2 for quarks) and S^^, stands for the spin 
part of the angular momentum operator. For a quark field ^{z) = q{z) 

S^^g = --cTf^uQ- (3.3) 
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In this paper we will be interested in the conformal transformations for the fields 'living' 
on the light-cone 

<!?{z) = <^{znf,) , n^ = 0. (3.4) 

One can check that the only remaining nontrivial generators are P+, D, M+_ and K_, 
where '+' and '— ' stand for the projection on n^ and on the alternative light-like vector 
n, fin = 1, respectively. We will further assume that the field $ is chosen to be an 
eigenstate of the spin operator E__|_, that is it has fixed spin projection s on the '+' 
direction^ 

S_+$ = is<l>. (3.5) 



L+ = Li -iL2 = -K_, 


L_ = Li + iL2 


Lo = ^(D + M_+), 


E = ^(D-M_+). 



For the leading-twist quark operators ( !2.2U| ), ( ^.21 ) s = -1-1/2 for each of the three quarks 
since S_+(^g) = +i{flq)/2. 

To bring the commutation relations ( p.l|) to the standard form it is convenient to 
consider the following linear combinations: 

-2P + , 

(3.6) 

The operators Lj form the so-called collinear subalgebra SL{2, M.) of the conformal alge- 
bra: 

[L+,L_]=2Lo, [Lo,L+]=L+, [Lo,L_] = -L_. (3.7) 

Most importantly, action of the group generators ( p.6|) on quantum fields $ (which is 
derived from (|3.2| ) by simple algebra) can be replaced by differential operators acting on 
the field coordinates and satisfying the same SL{2) commutation relations:^ 

[L_,<|.(^)] = -^$(^)^L_$(^), 

[L+, H^)] = [z^j-^ + 2j z^ ^{z) ^ L+$(z) , (3.8) 

A one-particle operator ^{z) is an eigenstate of the quadratic Casimir operator 

L2 = l2 + L^ + l2 = l2 - Lo + L+L_, [L2, LJ = 0, (3.9) 

[^{z)] = j(j - iMz), (3.10) 

•^ This property is automatically satisfied for leading twist operators which correspond to the max- 
imum spin projection; in the general case one should use suitable projection operators to separate 
different spin components, see e.g. p5[ |. 

^ Note that we use boldface letters for the generators acting on quantum fields to distinguish from 
the corresponding differential operators acting on the field coordinates. 
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where 

j = ^{l + s). (3.11) 

We will refer to j as conformal spin of $ in what follows. The remaining generator E 
counts the twist t = / — s of $: 

[E,<l>{z)] = ^-{l-s)<l>{z). (3.12) 

It commutes with all Lj and is not relevant for further discussion. 

It is helpful to have in mind that the operators La generate the projective (Mobius) 
transformations on the line in the '+' direction on the light-cone: 

, az + b 

z ^ z = , ad — be = 1 , 

cz + a 

<^{z) -^ $'(^) = {cz + rf)-2^$ (^^^) (3.13) 

with a, b, c, d real. The collinear conformal transformations of the three-quark opera- 
tors B{zi, Z2,Z3) defined in ( p.20|) and ( p.21| ) correspond to independent transformations 



( p.l3| ) for each of the fields; the group generators are given by the sum of one-particle 
generators acting on light-cone coordinates of the quarks: 

LaB{zi, Z2, Zs) = (Li,„ + L2,a + L3^a)B{Zi, Z2, Z3) , (3.14) 



where a = {0, +, — } and L^^a is the differential operator ( p.8[ ) acting on the argument 
of the A;— th quark Zk, k = {1, 2, 3}. For further use, we introduce two- and three-quark 
Casimir operators 

^ik — /_^ [Li^a + Lk^a) 5 i,k = 1,2,3, 
a=0,l,2 



L — 2_^ (-^1,0 + -^2,a + -^3,0) 
0=0,1,2 

= LI2 + Ll, + Ll - Ji(ji - 1) - J2{j2 - 1) - J3(j3 - 1). (3.15) 



The last three terms in the last line vanish for quark fields for which j = 1, see ( 3.11 
The two-particle Casimir operators L|^ can be written (for quarks) as 

L^ik = -didkZik, (3.16) 

where dk = d/dzk and Zik = Zi — Zk- Obviously [L^^, L^] = 0. 
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3.2 Brodsky-Lepage equations in the SL{2) covariant form 

The expected conformal invariance of the evolution equation for baryonic operators im- 
phes that the two-particle kernels Tiik commute with the generators of the SL{2) trans- 
formations La defined in ( p.l4|) and ( p.8|) . To show this, consider the following expression 
that generahzes both (lOSi) and (|]2^: 

Hi2B{zi,Z2,Z3) = Vauj{ai,a2)B{zi-aiZi2,Z2 + a2Zi2,zs), (3.17) 



where zu = zi — Z2 and the integration measure was defined in ( p.6|) . This operator 
has a simple meaning — acting on the three-particle nonlocal operator B{zi,Z2,Z3) it 
displaces the quarks with the coordinates zi and Z2 on the light-cone in the direction of 
each other. 

It easy to see that for this ansatz [Hu, L-] = [H12, Lq] = for an arbitrary function 
uj{ai,a2), whereas the condition [Hi2,L^] = leads to the following constraint: 

- — aiai + 2aiji]uj{ai,a2) = (^ — a2a2 + 2^2 J2 J ^("i, "2), (3.18) 



where a = 1 — a. Its general solution has the form 

dJ^ai, 0:2) = «i «2 V 



^n-2,2n-2,^ f ^ , (3.19) 

V"l"2/ 



with an arbitrary ip. However, remembering that the function y? should result from 
the calculation of the one-loop diagrams shown in Fig. |l| and must lead to nonsingular 
(bounded) operator H12, one may conclude that the form of ^ is almost uniquely fixed. 
Note that jk = 1 for all the three quark fields entering (p.20|) and (|2.21|) . Then, notice 



that the gluon exchange between quarks in Fig. |l|b amounts to the displacement of the 
two participating quarks along the light-cone and the function Lp must have a smooth 
behavior around ai = 0:2 = and ai = 0.2 = 0. These conditions leave us with the only 
choice p{x) = 1 and its substitution into ( p.l7| ) yields indeed the kernel ( p. 27] ). In a 
similar way, the 'vertex correction' in Fig. |l]a obviously corresponds to the displacement 
of just one of the quark operators and this leads to the second structure ip{x) = 6{x) 
which reproduces the two-particle kernel (|2.25|) . f\ 

Once conformal symmetry of the two-particle kernels is established, the group theory 
tells that Tiik may only depend on the corresponding two-particle Casimir operators 
L^^. To find the functional form of this dependence, one has to compare their action 
on a suitable basis of trial functions. The trick which we use below is general, and the 
calculation presents an example of the use of the 'dual basis' which is elaborated later 
in Sect. 3.5. 



^The second possible candidate (p{x) — S{1 — x) is ruled out since for ai + a2 — > 1 it gives rise to 
the operator B(zi ~ aizi2, zi — (1 — 02)^12, -^3) which becomes local in two quark fields. Such 'contact 
interaction' terms possess additional UV singularities and are not expected to appear. 
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For definiteness, let us find Hu as a function of Lfg- To this end, it is enough to 
compare their action on the homogeneous polynomials of two variables Zi and Z2: 

B{zi,Z2,Zs) > bn{Zi,Z2) 

which we choose to be eigenf unctions of the operator L^2 = —did2zf2- 

It is easy to see that the thus defined polynomials form an (infinite-dimensional) 
representation of the SL{2) group on which the operators L+ = Li^+ + L2,+ and L_ = 
Li^_ + L2- act as rising and lowering operators, respectively. It is thus sufficient to 
consider only the functions (polynomials) annihilated by L_ , or equivalently the highest 
weight of the representation, since all other eigenf unctions of L^2 can then be obtained 
by a repeated application of L^. Since L_ = —{di + 82) the latter condition is simply 
the translation invariance which leaves one with 

6(zi, Z2) = (^1 - Z2T = ^2, n = 0, 1, 2, . . . (3.20) 

An explicit calculation gives 

Ll2z'^2 = in + 2){n + l)z'^2, 

H^2^r2 = 2[^(n + 2) -^(2)1^^2, 

ni2z'l2 = ^/[in + 2){n + l)]z^2, (3-21) 

where iplx) = d\nr{x)/dx is the Euler ■i/'-function. To cast ( p.21| ) in an operator form, 
define J12 as a formal solution of the operator relation^ 

LI2 = Ji2{Ji2 - 1). (3.22) 

The eigenvalues of J12 equal ju = n + 2 and specify the possible values of the sum of 
two j = 1 conformal spins of quarks in the (12)-pair, cf. ( p.lO| ). Then 

H^2 = 2[^(Ji2)-V^(2)], 

^12 = l/[MJi2 - 1)] = I/LI2. (3.23) 



Substituting the representation (|3.23| ) into ( p.24|) and (|2.26|) one obtains the Schrodinger 



equation ( p.32|) for the three particles on the (light-cone) line with the coordinates zi, 
Z2 and Z3. The 'Hamiltonians' 7^3/2 and ?^i/2 entering this equation for different baryon 
states have a pairwise structure and are expressed in terms of the corresponding two- 
particle Casimir operators ( p.23| ). Furthermore, as we will elaborate in the next section, 
the Hamiltonian 7^3/2 possesses an additional 'hidden' symmetry: One can construct 



Since (3.22) is invariant under the substitution J12 — > 1 — J12 one has specify which one of the two 



formal solutions of (3.22) to choose; the simplest way to fix the solution is to take the one with larger 
eigenvalues. 
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an integral of motion (conserved charge) that commutes with 7^3/2 and with the SL{2) 
generators: 



Q = 2 [^12' ^23] = ^ ^1^2532:12223231 , 

[Q,L^] = [Q,n3/2] = 0. (3.24) 

Its presence makes the corresponding Schrodinger equation completely integrable and 
allows us to calculate the spectrum of the anomalous dimensions analytically by apply- 
ing a powerful technique of integrable models. The commutativity [Q,TC3/2] = is a 
consequence of the commutation relations between Q and two-particle Hamiltonians 

[^12, Q] = ^(^23 - ^31) , [^23, Q] = KLl - Ll,) , [HI,, Q] = i{Ll, - Ll,) . (3.25) 

The easiest way to prove these operator identities is to calculate both sides using the 
conformal basis of functions introduced below in Sect. 3.4 (see Eqs. ( p. 401) and (|3.41|) ). 



3.3 Conformal symmetry of the eigenfunctions 



Equations ( p. 231) define the Brodsky-Lepage evolution kernels in the most general form, 



independent on the representation. The particular choice of the SL{2) generators ( |3.8| ) 
corresponds to the evolution of the nonlocal operator B{zi, z^^ 2:3). As we have argued in 
Sect. 2.2, diagonalization of the evolution equation for baryon distribution amplitudes 
rather involves solution of the corresponding Schrodinger equation for the local operators, 
or, equivalently, their coefficient functions. This corresponds, formally, to going over to 
a different representation, and it is important to realize that the action of the generators 
of the coUinear conformal group on the elementary fields and on the coefficient functions 
of local operators defined through Eq. (|2.31|) is not the same. By requiring 

[2±,0^(5l, ^2, 93)]5(2i, Z2, 23) = ^(9i, ^2, 93))[L^,o5(zi, ^2, 23)] (3.26) 



2i=0 



2i=0 



one finds the following 'adjoint' representation of the generators acting on the space of 
coefficient functions '^{xx.Xi-.x-i): 

^fc,o^(a:i) = (xfcSfc + 1) ^{xi) , 

Lfc,+^(xi) = -Xfc^(xi) , 

Lfc,_^(xi) = [x^dl + 29fc) ^{xi), (3.27) 



where, in order to maintain the same commutation relations ( |3.7| ), we have defined L_ 
as the adjoint to L_|_, and vice versa. To simplify the notations, in what follows we drop 
the 'hat' from the generators in the adjoint representation, which, hopefully, will not 
yield confusion. Thus, the two-particle Hamiltonians entering the Schrodinger equation 
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for the coefficient functions, ( p.32| ), are given by the same operator expressions ( p.23| ) 



but with the SL{2) generators defined by ( p.27| ) 



As usual in quantum mechanics, symmetry of the Hamihonian imphes that the eigen- 
states are degenerate: applying the SL{2) generators to a particular eigenstate '^N,q one 
arrives at a yet another eigenstate with the same value of energy E^^q- It is then natural 
to parameterize the eigenstates '^N,q by a complete set of mutually commuting conserved 
charges. The conformal symmetry allows to identify two such quantum numbers: the 
total conformal spin, L^, and its projection, Lq, which are common to both Hamiltonians 
7^1/2 and 7^3/2- 

The construction of conformal eigenstates is fully analogous to the construction of 
the eigenstates of angular momentum in standard textbooks on quantum mechanics, 
with the 0(3) symmetry replaced by SL{2). We require that ^(xj) should diagonalize 
simultaneously two integrals of motion 

L^m{x^)=h{h-l)^!{xi), Lo^(x,) = (Ar + 3)^(x,). (3.28) 

Here, the first condition defines the conformal spin of the state, /i, and the second one 
follows trivially from the fact that \l/(xj) is a homogeneous polynomial of degree A^ 
in three variables Xj. Assuming that there exists a positive definite scalar product on 
the space of the coefficient functions (see (Eq. (|3.31|) below) one can easily prove that 
eigenvalues of L+L_ are non-positive. From the definition L^ = Lq{Lq — 1) + L+L^ it 
then follows that h < A^ + 3. Moreover, the eigenstate with the largest conformal spin 
h = N + ?) has to be annihilated by the lowering operator L_[] 

L_ ^(°) (xi) = , h = N + 'i (3.29) 

and is, thus, the highest weight of the representation. All other states can be obtained 
from the highest weight by a repeated application of the rising operator L_|_ which acts 
trivially on the coefficient function 

^(")(x,) = L!^^W(a;,) = (-l)"(xi + X2 + x^Y^^''\xi) , 

L^¥'^\x,) = {h + n)¥''\xi) (3.30) 

and amounts to 'dressing' the corresponding local operator by the n— th power of a total 
derivative. These states form an infinite dimensional representation of the SL{2) group 
of a positive discrete series labeled by the integer conformal spin h = N + 3. They all 
have the same energy and, being substituted into (|2.19| ) and (|2.35| ), lead to identical 



contributions to the baryon distribution amplitude due to the condition xi + X2 + X3 = 1. 
For this reason we can neglect such states altogether and impose ( |3.29| ) as an additional 
constraint on the solutions of the Schrodinger equation ( |2.32| )p|. 

^ L_|_ and L_ act as the rising and the lowering operators in the (infinite dimensional) representation 
labeled by the spin h, respectively, so that if Lq'^ = /iq^ then LoL±'^ = {ho ± 1)^. 

^As familiar from quantum mechanics, the highest weight states exhibit additional symmetry. In our 
case, it is easy to find that the local operator corresponding to the highest weight transforms under the 



SL{2) transformation according to (3.13) 
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Note that the conformal spin h of the three-quark state which satisfies the highest 
weight condition ( |3.29|) is related to the total number of derivatives. As a consequence, 
conformal operators with different A^ do not mix with each other under renormalization 
since they belong to different representations of the collinear conformal group. This 
condition is yet not sufficient to diagonalize the evolution equation since, as will become 
clear in the next section, for fixed N there exist A^ + 1 different conformal operators 
mixing between which is allowed by conformal symmetry and exists, in general. The size 
of the mixing matrix is, however, reduced from A^(A^ + l)/2 to A^ + 1. The impact of 
conformal symmetry is that one can eliminate all mixing with operators containing total 
derivatives. 

It is straightforward to check that the SL{2) generators La as well as the Hamiltonians 
in ( |2.26| ) and ( ^.24 ) are Hermitian with respect to the 5'L(2)-invariant scalar product: 



(^l|^2) = 120 VxXiX2X3%iXi)^2ix,). (3.31) 

Hermiticity implies that the scalar product of two eigenfunctions with different eigen- 
values vanishes, i.e. coefficient functions of any two operators that do not mix under 
renormalization are orthogonal with this weight function, cf. ( 2.34| ). 



3.4 The conformal basis 

To find the general solution of the evolution equation ( ^.32| ) with the Hamiltonian 7i 
given in ( p. 241 ) or ( p.26| ) it proves convenient to decompose the eigenfunctions \E'(xj) over 



,{12)3, 



a suitable basis of functions Wn (a^i) having the same conformal properties as ^{xi): 

N 

^(xi,X2,a;3) =5^^"M„/-i^ii')3(xi,X2,X3). (3.32) 

n=0 

Here, the factor z" is inserted in order that the coefficients u„ are real, as will become 
clear in the next section. The numerical factor 

= ^^f^iylM/^ - 1) - Ji2(ji2 - 1)] (3.33) 

is included for later convenience and the notations \l'n and ji2 = n+2 will be explained 
below. 

Aim of this section is to construct such a basis. To this end, we require that the 
functions vl'^^^)^ have the same conformal properties as \E', that is 

Lo *r)'(x,) = h¥f^^x,) , h = N + 3, 

L' vl/a2)3(a;^) = h{h - 1)^1/^)3(3..). (3.34) 
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The second-order differential equations Eq. ( p.34| ) do not specify the set of polynomi- 
als ^n uniquely, but rather allow to choose them as a linear combination of (A^ + 1) 
solutions with arbitrary coefficients. To fix these coefficients one has to supplement 
Eq. (|3.34| ) by some additional condition. The traditional choice [§ is to expand \E'(xj) 
over the set of Appell polynomials ^6) (for xi + X2 + 0:3 = 1): 

An,N-n{Xl, X2) ~ [XiX2(l - Xi - X2)Y^ dld^-"" x\^''x\^''-''{l - Xi - ^2)^+^. (3.35) 

In this way, solving the evolution equation ( p.32| ), one is left with a complicated (A^ + 
1) X (A^ + 1) mixing matrix for the coefficients in front of Appell polynomials with the 
same A^ but different n, which does not have any obvious structure. This basis is also 
inconvenient for calculations since Appell polynomials with different values of n are not 
mutually orthogonal. 

The expansion in Appell polynomials is, however, not warranted and in this paper we 
suggest a different basis which is orthonormal and better suited for the solution of the 
evolution equation. To this end, we require that in addition to ( 3.34|) the polynomials 



^n {xi) (n = 0,...,A^) should diagonalize the two-particle Casimir operator in the 
channel defined by the (12)-quark pair: 



Ll,^>^^^^\x.)=n2{jl2-l)^^n^^\x 



*7? 



J12 = n + 2 , Q<n<N . (3.36) 

The particular choice of a quark pair is of course arbitrary and we might use, e.g., L\^ 
for the same purpose. In this way one obtains a different basis of functions '^n \Xi) 
that are linear related to Wn (xj) through the Racah 6j— symbols of the SL{2) group 
(see Appendix A): 

N 

^T\^^) = Y.^kn^n'^\^i) ■ (3-37) 

n=0 

Here, the superscript indicates the order in which the tensor product of three SL{2) 
representations has been decomposed into the irreducible components. 

The solution of the combined Eqs. ( p. 34 ) and ( 3.36| ) can be obtained either solving 



the corresponding second-order differential equations explicitly or making use of the 
conformal OPE. The result reads (in a certain convenient normalization): 

*r(-.) = (/v+»+4)(x.+..r(x.+x.+.3)"-pS"'( "; : "■ : T Wi"-^) 

\Xi + X2 + X3 J \Xi-\-X2/ 

(3.38) 
where Cj (x) and P^ (x) are Gegenbauer and Jacobi polynomials |^, respectively. 



Note that each function ^]y ^ is specified by a pair of integers A^, n which are related 
in a obvious way to the total conformal spin of the three-quark operator h = N + 3 and 
the conformal spin of the (12)-pair J12 = n + 2, respectively. In what follows we often 
drop the subscript W if it is clear from the context. 
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3.4.1 Properties of the conformal basis 

The foUowing features of the new basis are especially important. 

First, the functions ^jvn (^0 ^^^ mutually orthogonal with respect to the SL{2) 
scalar product ( p.31|) 



^^Sf I^IS) = 120 IVxx,x,x,^>fhx.)^^^^\ 



Xi 



SmnS. 



MNOmn 



QOfn 

2N + 5' 



(3.39) 



The integration measure JVx is defined in ( p.6|) , /7v,n is given in ( p.33|) and the factor 
120 is introduced in order that fVx ■ 120xia;2a;3 = 1. 

Second, action of the Casimir operators of the coUinear conformal group in this 
basis is rather simple. By construction, \l/jv„ (xj) diagonalize L^ and Lfg whereas the 
remaining two two-particle Casimir operators turn out to be three-diagonal: 



^L<")^( 



X'i 



Ll,<-)3(x. 



In 
In 



,(12)3 



(n + 1) 
1 

(n + l) 



CT(^0 + 



eT(x.)+ 



2n + 3 



(n + 2)(n + l) 

2ra + 3 
(n + 2)(n + l) 






(n + 2) ^"+^ ^"'^ 



In + 2) ^"+^ ^"^^ 



(3.40) 



This property turns out to be crucial for simplification of the evolution equation. In 
particular, using the definition ( p.24| ) one finds that the operator Q can be represented 
in the conformal basis by a (A^-|- 1) x (A^-|- 1) matrix with only two subleading diagonals 
nonzero 



Qv&r^^i 



Xi 



ifn 



^n+1 I"''* 



't'Si^, 



(3.41) 



Finally, the factorized form of \l/„ (xj) as a product of polynomials depending sep- 
arately on s = (xi — X2)/(xi + X2) and t = (xi + X2 — X3)/(xi + X2 + X3) is conve- 
nient for applications. Note that the integration measure is also factorized: jT>x = 
/_>(t + l)/4/_\rf.. 

3.4.2 Special cases 



The definition in ( |3.38|) is valid for arbitrary xi,X2,X3. One important special case is 
Xi + X2 + X3 = 1 which corresponds to the expansion of distribution amplitudes so that 
Xi = {xi,X2,X3} can be identified with the set of quark momentum fractions. For this 
case we obtain a complete set of polynomials 



*!if (^.) = {N + n + A)ix, + x.rPt:'''^ (2X3 - 1) Cf ^ 



Xi - X2 

Xi + X2 



(3.42) 



which, as we are going to argue, are much superior for studies of the three-particle 
distribution amplitudes as compared to Appell polynomials. 
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Another important case is xi + X2 + xs = which corresponds to neglecting contribu- 
tions of all operators containing total derivatives. This choice is relevant if, for example, 
one considers only forward matrix elements. It also allows to abstract from unnecessary 
'kinematical' complications related to the conformal symmetry and consider the dynam- 
ical mixing problem in the most pure form. Note that the basis functions (p.38|) become 
very simple: 



^1^^)=^^^ 



= w{N, n) (xi + x,rcT -^-^ , (3.43) 



Xi - X2 



where 



(2iV + 4)! _ , ,,.-,,, V{2h-l) 



wiN,n) = i-ir~ \^, , : , , „ : „„ = -(-1)^ 



{N + n + 3)\iN-n)\ ' ' r(/i + ji2 - l)r(/i - J12) ' 

(3.44) 

Instead of the full coefficient function \E'(xj) one can consider the function of one variable 

'^{x) defined as 



X\X2X-i'^{Xi 



-(xi + X2)^+3vl/(^^^), (3.45) 



,(12)3, 



SO that if \E' is expanded in the basis of \E'n (xj) with the coefficients as in Eq. (|3.32| ), 
then 

^{x) = ^^Y.V^u^f-^w{N,n)CT{x). (3.46) 

n=0 

Note that although \E'(x) was obtained from the coefficient function \I'(xj) by reduction 
to the subspace xi + X2 + X3 = 0, it contains all nontrivial dynamics of the problem. If 
^(x) is known, then the full function \&(xj) of three variables can easily be recovered 
through its expansion ( p.32| ) since 



„ _,n, 1^;v(^-^+l)!(iV + ^ + 4)! /-^ 



(2Ar + 4)! 



I dx^{x)C^J'^{x). (3.47) 



In physical terms, existence of such a relation is a consequence of the triangular struc- 
ture of the mixing matrix with the operators containing total derivatives, familiar from 
studies of meson distribution amplitudes. Similar to the latter case, it is sufficient for cal- 
culation of the anomalous dimensions to consider forward matrix elements of three-quark 
operators for free quarks. After this is done, the coefficient functions of multiplicatively 
renormalizable operators can be obtained from ( |3.47| ), ( |3.32| ). 



The algebraic structure of this connection is, however, complicated, which can be 
traced to the fact that the lowering operator L_ is nontrivial in the 'adjoint' represen- 
tation ( p.27|) . As a consequence, there exists no simple way to resolve the constraints 



imposed on the form of the function \I'(xj) by the highest weight condition L_\I' = 0. 

In what follows we suggest an alternative basis in which the eigenf unctions have a 
much simpler form that is useful in some applications. 
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3.5 The dual conformal basis 

Once the evolution 'Hamiltonians' 7i are written in terms of the SL{2) generators, one 
can abstract from the 'physical' Hilbert space spanned by the coefficient functions '^{xi) 
and try to find an equivalent representation of the SL{2) group with simpler properties of 
the highest weights. The calculation made in Sect. 3.2 suggests that the conformal sym- 
metry properties of polynomials $(^j) of the light-cone coordinates Zi might be simpler 
than polynomials \l/(xi) of the momentum fractions since according to ( p.8| ) the highest 
weight condition L_^{zi) = translates to the translation invariance of ^(zj)^. The 
translation invariance, combined with the restriction to homogeneous polynomials ^{zi) 
of degree A^ in three light-cone coordinates Zi, implies that ^{zi) essentially reduces to 
a polynomial of degree A^ of a single variable, times a simple overall factor: 

\ ^ ^{z,) = {z, - z,r ^ r-^^] . (3.48) 

L_$(2,) = J \Z1-Z2J 

Note similarity to, and at the same time difference with Eq. (|3.45|) defining the coefficient 
function of one variable \l/(x) for the special choice of momentum fractions: In both 
representations the conformal symmetry allows one to reduce the evolution equation 
involving three variables to an equation involving a function of a single variable — \1/ (x) or 
$(2;), respectively. [^ At the same time, while the one- and the three- variable descriptions 
are essentially equivalent in position space thanks to the translation invariance, the 
relation between \l/(x) and \l/(xj) appears to be much less transparent, see Sect. 3.4. 

The easiest way to construct the basis of polynomials in position space explicitly is to 
identify them with suitable correlation functions in a certain two-dimensional conformal 
field theory. Let Ox^, be a local conformal operator with spin h = N + ?> corresponding to 
the coefficient function \l/(xj) so that it is transformed as an elementary field with spin 
h under the projective transformations ( |3.13|) . In so far as only these transformation 
properties are important, we can replace formally the quarks by free scalar fields 0(^fe) 
with the same conformal spin jk = 1: (^^(0 = ^(^i; ^^2, <93)0(^i)0(^2)0(^3)|t.=t- In the 
terminology of conformal field theories such operators are called quasiprimary fields. 
Correlation functions of them with elementary fields are known to satisfy the conformal 
Ward identities which take the form of the highest weight conditions, (p.28|) and (|3.29| ), 



that we are looking for. This suggests to define the polynomial ^{zi) as dual to the 

^ The basis of the functions ^{zi) is dual to the conformal basis ^(xi) in the same sense as the light- 
cone coordinates Zi are dual to the light-cone momentum fractions Xi. The 'physical' coordinate space 
distribution amplitudes were introduced in Rcf. ||2^: these are states which diagonalize the lowering 
operator L-VP = —ip^ and thus resemble coherent states in the standard field theory terminology. 
Here p has a physical meaning of the momentum of the hadronic 'wave packet' propagating along the 
light-cone. 

-•^^As we will show in Sect. 4.3, these two functions are related to each other through the duality 
transformation. 
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coefficient function \l/(xj) by the following correlation function: 

$(Z1,Z2,^3) = wlwlwl{0\O<i,{0)<f){wi)<P{w2)(l){w3)\0) 



Zi = l/Wi 



vi/(%,%,%)n(i-^fc6) 



fc=i 



€fe=0 



(3.49) 



where we used the expression for a propagator of free field (O|0(u7)0(O)|O) = w^'^. By 
construction, the SL{2) generators have the standard representation (|3.8|) on the space 
of dual polynomials and it is straightforward to verify that $(-2j) defined in this way 
satisfies the conditions ( p.48|) . 

The two polynomials $(2:1, 2:2, Zs) and '^{xi,X2,x^) have the same degree N and are 
related to each other by the Mellin transformation [^] 



^{zi,Z2,z^)= I Wdtktke *'=^(ziti, 2:2^2, 2:3^3) 



that amounts to the redefinition of the coefficients c 



ni,n2,n[i ^n-i,n2,nz 

„ ~ni -"2 ^nz 

n2,nz ^ni,n2,n3^1 ^2 ^3 • 



(3.50) 

>i + l)!(ri2 + 



l)!(ra3 + 1)! in the polynomial '<^ {zi, Z2, z^) = ^^^ 

The SL{2) invariant scalar product on the space of coefficient functions (|3.31 
equivalently be rewritten as 

r(2iv + 6),,,. 



can 



(^l|*2) 



r(6) 

r(2A^ + 6) 

r(6) 



'^l{d,,,d,^,d,,)^2{zi, Z2, Z3) 
'^2{dz^,d,^, d,,)<^l{zi, Z2, zs) 



Zk=0 



Zk=0 



(3.51) 



where $1(2) is a dual of \l'i(2). 

The two representations for the eigenfunctions, ^{zi) and \l/(xj), are equivalent from 
the point of view of diagonalization of the evolution equation: they give rise to the 
same energy spectrum and are related to each other through the transformation (|3.5CI| ). 
However, the use of the dual representation can be advantageous due to the particular 
simple structure ( |3.48| ). 

Applying the transformation ( |3.5CI| ) to the both sides of ( ^.321 ) one can construct 
the dual conformal basis ^'^^'^^'^ (zi) . The functions $7Vn ^^^ be defined as translation- 
invariant homogeneous polynomials of three variables which diagonalize the Casimir 



operator LI2 (i^i the standard representation (|3.8| )): 



,(12)3 



^N,n i^i) = (^1 - ^2) ^N,n{^) ' ^ 

^?2<''^'(^.) = (^ + 2) (n + 1X12)3 (^. 



^3-^2 

Zl- Z2'' 



(3.52) 
(3.53) 



with n = 0,...,N. Solving the last condition one gets the explicit expression for the 
functions {pN,n{z): 



<PnA^) = fr 



{N + n + 3)\{n + 1)\ 
(2n + 2)! 



^N-n 



Fi 



n- N,n + 2 
2n + A 



-1 



(3.54) 
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which defines <fN,n{z) as a polynomial of degree A^ — n in 2;. Here the normalization 
is such that the polynomials ( |3.52|) and ( |3.38|) are related to each other by the Mellin 



transformation (|3.5CI|) . The decomposition of the dual eigenfunction $(-2j) over the dual 
basis has again the form (|3.32|) with the same coefficients m„ 



N 



^Zi, Z2, Z3) = {Zi - Z2)^ Yl ^"^"/„ V7V,n(2;), (3.55) 



n=0 



with z defined in ( |3.52| ). 

It is clear that the linear algebraic relations ( |3.37| ), ( |3.40| ) and ( p.41| ) satisfied by the 
polynomials \E'*^^^)^(xj) remain valid for the dual polynomials $^^^'*^(zj) provided that one 
changes the adjoint representation of the SL{2) generators (|3.27|) to the standard one in 



The function ipN,n{z) has two indices corresponding to the total conformal spin of 
the system h = N + 3 and the conformal spin ji2 = n + 2 in the subchannel (12). In the 
sequel we will need the asymptotic behavior of this function in the limit when any two 
of the coordinates Zi coincide, or equivalently z = 0, 1 and oo: 



^_„(iV + ^ + 3)!(n + l)! 



^nA^) = fnz ^2n + 2)! 

^nA^) '= /n(2n + 3)(iV + l)!, (3.56) 

^nM '= /„(-l)^-"(2n + 3)(iV+l)!, 
where the leading terms are kept only. 

4 Integrability 



As was explained in Sect. 3.2, the Brodsky-Lepage evolution equations ( p. 32 ) have the 



form of the Schrodinger equations describing a three-particle system with three degrees of 
freedom which we can choose as quark momentum fractions Xi or coordinates Zi depend- 
ing on whether the 'physical' or 'dual' representation is usecQ for the 'wave functions'. 
Either way, the scale dependence of baryon distribution amplitudes in QCD corresponds 
to a one-dimensional quantum mechanical 3-body problem with very peculiar Hamil- 
tonians, ( p.24| ), ( p.26| ) and ( p.23| ), determined by the underlying QCD dynamics. The 



conformal symmetry allows to trade two degrees of freedom for two quantum numbers 
corresponding to the total conformal spin LF' and its projection Lq after which one is left 
with one degree of freedom described by either the set of coefficients Un in the conformal 
basis ( |3.32D or, equivalently, a function of a single variable ( |3.45| ) or ( ^.48] ). The orig- 
inal 3-body Schrodinger equation is reduced, accordingly, to a (comphcated) one-body 
problem which is in general not possible to solve analytically for arbitrary A^. 



^^Note that the expressions for 7i", 'W in Eqs. (2.25), (2.27) correspond to the dual representation. 
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The crucial observation is that the Haniiltonian 7Y3/2 (but not 7Y1/2) proves to be 
completely integrable: The operator Q defined in ( |3.24|) commutes both with the Hamil- 



tonian and with generators of the SL{2) group. The eigenvalues of Q thus provide us 
with the third quantum number allowing to specify completely the three quark states 
with maximal helicity. Existence of the nontrivial 'conserved charge' implies that 7^3/2 
is a (complicated) function of two and only two mutually commuting operators Q and 
L^. Therefore, instead of solving the Schrodinger equation ( |2.32| ) directly, one can solve 
much simpler equations ( |3.28|) supplemented by the additional condition 



Q <!>{zi) = -id,,d,^d,^zi2Z23Z3i'^{zi) = q $(zi), (4.1) 

Q'^ixi) = i {d^, - (9^.J (9^2 - 9^3) (9^.3 - (9^. J XiX2X3-^{xi) =q'^{xi) (4.2) 

in the 'dual' and the 'physical' representations, respectively, and find the spectrum of 
the Hamiltonian 7^3/2 = 7i3/2(-^^,Q) by replacing the operators by their corresponding 
eigenvalues. 

Remarkably enough the Hamiltonian 7^3/2 is well known from integrable generaliza- 



tions of the Heisenberg spin magnet models |jT8[. Indeed, an inspection shows that the 
SL{2) generators (p.8|) for quarks with conformal spin jk = 1 can be interpreted as 
Lorentz spin s = — 1 operators. In this way, we may consider the Hamiltonian 7^3/2 as 
describing the system of three interacting spins each acting on its internal space labeled 
by the coordinates Zk- These spins carry the index of the corresponding particles and 
form a one- dimensional spin chain with three sites. This system coincides identically 
with the celebrated one-dimensional XXX Heisenberg spin magnet of noncompact spin 
s = — 1 for which powerful Quantum Inverse Scattering methods have been developed 
and a lot of results are available [^-@- Aim of this section is to elaborate on this 
connection and adapt the existing results to the present context. Some new results will 
be presented as well. 

4.1 The master recurrence relation 

By construction of the conformal basis, the eigenfunctions ( |3.32| ) and ( |3.55| ) obey the 
conditions ( |3.28| ) for arbitrary coefficients Un- Using Eq. ( |3.41| ) it is easy to derive that 
the equation Q\l/ = g\l/ is equivalent to the following three terms recurrence relation for 
the coefficients u„, {n = 0, . . . ,N): 

qUn = fniUn+l + Un-l) , (4.3) 

with the 'boundary' conditions 

u_i = un+1 = 0, (4.4) 



which follow from the properties of the coefficients /„ (|3.33| ). The overall normalization 
of Un is arbitrary and we choose for simplicity 

Uo = l. (4.5) 
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The recurrence relations ( [4 .31 ) represent the system of A^ + 1 hnear homogeneous 
equations on the coefficients u^- Solution of this system is equivalent to diagonalization 
of a (A^+1) X (A^+1) matrix with only two subleading diagonals nonzero. The consistency 
condition for this system translates to the characteristic polynomial of degree iV + 1 in 
q whose zeros define the A^ + 1 quantized values of q. 

It follows from the recurrence relations ( |4.3| ) that Un{q) (n = 0, ..., A^) form a system 
of (semiclassical) orthogonal polynomials in a discrete variable q. Then, the boundary 
condition uj\f+i{q) = implies that A^ + 1 quantized values of q have the properties of 
roots of orthogonal polynomials, that is, they are real and simple, for different A^ the set 
of quantized q are interlaced. The completeness and orthogonality conditions for this 
system are given by the Cristoffel-Darboux relations PB 



y" 7 rTUn{q)Un{q) = ^qq' , (4.6) 

— jnUJyq) 



n=0 



y^T — T^Un{q)um{q) = S^m, (4.7) 

g fnUj{q) 

where u{q) = u^{q)dqU]^^i{q) and in the second line the summation goes over A^ + 1 
quantized q. 

The orthogonal polynomials Un{q) have an obvious parity property 

Un{-q) = i-lTuniq), (4.8) 

where from it follows that all nonzero eigenvalues of q come in pairs: If q is an eigenvalue, 
then —q is also an eigenvalue, U]sr^i{—q) = 0. In addition, for any even N there is a 
single eigenvalue g = and the corresponding coefficients are given by 

U2k{q = 0) = (-1)'% , U2k-i{q = 0) = (4.9) 

for k = l,...,N/2. 

4.2 Permutation symmetry 

The Hamiltonian 7^3/2 is explicitly invariant under cyclic permutations of the three 
particles. We define the generator of the corresponding discrete transformations V as 

V $(^1, Z2, zs) = $(^2, zs, z,) , [V, 7^3/2] = [V, L^] = [V,Q] = 0, (4.10) 

where, for definiteness, we have chosen to use the dual (coordinate space) representation. 
Because of the symmetry, the eigenfunctions of H3/2 can simultaneously be chosen as 
eigenstates of V: 

V^Zi,Z2,Z3) = ^Z2,Z3,Zi) =e^Zl,Z2,Z3) (4.11) 
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with 9 = 6{N, q) being a function of quantum numbers. Since P^ = 1, the possible 
eigenvalues 6 are given by three different cubic roots of unity: 

9 = e-^^(^-) , <p = 0,f,f. (4.12) 

In addition, 7^3/2 is symmetric under permutations of quarks in the (12) pair 

Vu ^z^, Z2, zs) = ^Z2, z^, z^) , [Pi2, ^3/2] = [^2, L'] = . (4.13) 

This implies that the eigenstates can be chosen to possess a definite 'parity' ^12 = ±1. 
In fact, the spin and isospin symmetry of the physical baryon distribution amplitudes 
introduced in Sect. 2 lead to their definite parity properties, see Eqs. ( p.3| ), ( p.l5| ), so 
that expansion in parity eigenstates is natural. 

One should stress that the operators V and 7^12 do not commute and therefore the 
eigenstates of V do not have, in general, definite parity, and vice versa. Nevertheless, 
the symmetry of the Hamiltonian under both V and P12 immediately implies that the 
eigenvalues of if 3/2 with 6^1 have to be (at least) double degenerate.^ 

Integrability of 7^3/2 alias existence of the conserved charge Q increases the symmetry, 
so that the 9 = 1 eigenstates turn out to be double degenerate as well, apart from the 
singular state corresponding to g = 0. To show this, note that P12 anticommutes with 
Q: 

V12Q = -QVi2- (4.14) 

Since the Hamiltonian 7^3/2 commutes simultaneously with Vu and Q, it should be an 
even function of Q and therefore the levels corresponding to nonzero q and —q have the 
same energy and are double degenerate. 

It follows from ( ^.14] ) that permutation of quarks transforms an eigenfunction of Q 



into another eigenfunction with the opposite value of q and the same value of the energy: 

$,(^2, Zi, Z3) = ^-q{Zi, Z2, Zs) = ^liZi, Z2, Z3) , 

^g(x2, xi, X3) = ^-^(xi, X2, X3) = ^*(a;i, X2, X3) . (4.15) 

These relations are an obvious consequence of ([4.1|) and ( |4.2|) . For the corresponding 
functions of one variable one gets: 

Re$(l-^) =Re$(2), Im$(l - ^) = -Im$(z) , (4.16) 

Re ^(-x) = Re ^(x) , Im ^(-x) = -Im ^(x) . (4.17) 



^^ To show this, consider an eigenstate of 7^3/2 which is simultaneously an eigenstate ofV: 7^3/2$ = 
£$ , V^ = 9^ . Acting on the first equation by 7^12, one gets H3/2'Pi2^ = £7^12$ , so that either 7^12$ is 
an independent eigenstate with the same energy, or it is proportional to $: 7^12$ = p^ and is therefore 
a parity eigenstate with p = ±1. In the latter case, applying the identity V12PP12 = T'^ to <i> one gets 
9 = 9^ where from necessarily 9 = 1. 
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Eqs. ( [4.15| ) suggest that real and imaginary parts of the complex eigenf unctions $(-2i), 
"^{zi) have definite parity with respect to the 7^12 permutations. Define 

^^{x,) = ¥+\x,)+t^'f\x,) (4.18) 



with the real functions 



^i-'\x,) = ^[<il,ix,) + ^.,ix,)], (4.19) 



<i'i-\x,) = -'-[^,ix.d-^-,ix,)]. (4.20) 

Then \1'+ (^^) is even (odd) with respect to permutations of the two first arguments: 

¥^^\xi,X2,X,)=±¥^^\x2,Xi,X3). (4.21) 

We recall that the eigenstates "^q correspond to the same value of the energy but, in 
contrast to \E'±g, they do not correspond, in general, to any definite eigenvalue 9{N,q). 
The eigenvalues ^(A^, q) of the cyclic permutation operator V can be expressed in 
terms of the solutions of the recurrence relation. To this end, substitute ^{zi) in ( |4.11| ) by 
its expansion in Eq. ( p.55| ) and take into account that the cyclic permutations correspond 
to the following transformation rules for the coordinate ratio z = {z^ — Z2) / {zi — Z2): 

zAl--A Uz. 4.22 

z 1 — z 

This gives 

TV . . N 

^-l)N^Nj2i-uJ-'^Jl--\ = J]r«„/-Vn(^), (4.23) 

n=0 ^ ^ n=0 

which has to be valid for an arbitrary real z. Consider the limit 2; — > cxd or, equivalently, 
zi — Z2 ^ 0. Taking into account Eqs. ( p.56[ ) we compare the leading asymptotics of the 



both sides of (|]23]) to get 

•'" (iV + 2)(iV + 3)^ "''" ' Eloi"«»(«)(2n. + 3) ' 

where the second equality follows from the identity 9 = 9* /9 and reality of the coefficients 
Un as defined by the recurrence relation with real coefficients. Comparing ( [4.24| ) with 
( [4.12| ) we end up with 

" N 

^2"M„(g)(2n + 3) . (4.25) 



0(A^, g) = 2 arg 



.n=0 

In terminology of integrable models this expression defines the quasimomentum corre- 
sponding to the wave function ^(zi). 

Since 0(A^, q) takes a discrete set of values ( |4.12|) , Eq. ( |4.25| ) suggests that eigenvalues 



of Q can be parameterized by an integer number (, and belong to a one-parametric 
family of curves, q = q{N,i). We will elaborate on the physical interpretation of such 
trajectories in Sect. 4.4.3 and construct them explicitly in Sect. 4.5 using the WKB 
expansion. 
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4.3 Duality 

The Hainiltonian 7^3/2 possesses a duality symmetry |31|| which allows to establish the 
equivalence between coefficient functions ^(xj) at a;i+a;2+a;3 = and the dual coefficient 
functions $(-2i). We recall that both functions are related to each other through the 
integral transformation ( p. 501 ) that maps the momentum fractions Xi into the light-cone 
coordinates Zi. 

As a hint, observe that changing the variables in ( [4.2| ) as xi — > Zi2, x^ — ^ -^23 and 
2^3 ~^ -231 one can formally cast it into the form of ( [4.1| ). In order to establish a formal 
equivalence, define the duality transformation S as 

^k *■ oXkO = Xk 3^fc+l) 

d., - 5.,,, - S {d,, - 9.,,J S-' = -9.,,,, (4.26) 

'i'{xi,X2,Xs) -> S'^{xi,X2,X3) = ^(xi2,a;23,a;3i), 

with Xfc+3 = Xk and k = 1,2, 3. Here, the second relation follows from the remaining 
two. It is easy to see that the constraint xi + X2 + 0:3 = is required as the consistency 
condition for these transformations. Using the definition and taking into account that 
'if{xi,X2,X3) is a homogeneous polynomial of degree N it is easy to check that 

S^XkS-^ = -3xk+i, (4.27) 

S^ ^(xi, X2, X3) = ^(-3x2, -3x3, -3xi) = (-3)^^(x2, X3, xi) , (4.28) 

which allows to express S'^ in terms of the cyclic permutation operator, V, and the SL{2) 
generator Lq as 

S^ = {-3)^"-^V . (4.29) 

Thus, the operator of the duality transformation S is formally proportional to the square 
root of the operator of cyclic permutations. 

Applying the transformation ( [4.26|) to the conserved charge Q in the adjoint repre- 



sentation (^4.2|) we find the expression 

SQS~^ = -id^^d^^d^^xi2X23X3i = Q , (4.30) 

which, after the replacement of the momentum fractions by the coordinates, x^ -^ z^, 
coincides with the operator Q acting on the dual coefficient functions, ( |4.1| ). For clarity, 
we have restored the 'hat' to indicate the adjoint representation, cf. ( ^.27| ). In the similar 
way, one can check that on the subspace xi + X2 + X3 = 

SL^S-^ = S 5^(xfc9,, + l)S-' = ^(x,9,, + 1) = Lo , (4.31) 

k k 

SL^S-' = -S Y,id.. - d.,fxkX,S-' = -Y,d.,d,^{xk-x,f = L'. (4.32) 

j>k j>k 

29 



In other words, the duahty transformation S maps the conserved charge Q and the 
SL{2) generators, Lq and L^, in the standard and the adjoint representations, ( p.8|) and 
( p.27| ), one into another. Since in both descriptions they form a complete set of mutually 
commuting operators, it follows that the eigenfunctions must transform one into another 
as well, up to a numerical factor: 



= C- <^N,g{zi,Z2,Z3) (4.33) 



"^ N,q{xi, X2, X3 

with C = C{N, q) being a normalization constant. Its value can be found by examining 



the asymptotics of the both sides as zi — zs —> 0, or equivalently z ^ 1. Using ( |3.32| ) 
( p.43| ) and (|3.55|) it is straightforward to get 



^■1 ,^9uo {N+iy.{N + 2)\iN + 3)\ 

^ ~ 2uN (2iV + 3)! • ^ '^ ^ 

The duality relation (|4.33|) is highly nontrivial and it is easy to see that this relation 
does not hold for the basis functions \l/(^^)^ and $(^^)^. The reason for this is that the 
defining relations (|3.36| ) and (|3.52|) are not mapped one into another by the duality 



transformation since 

SLl^ S'^ = -S {d^^ - d^^fxiX2 S'^ = -dlxi2X2i, ^ -d^^d^^xj^ . (4.35) 



This also explains why the expansions ( |3.46| ) and ( |3.55| ) involve the same coefficients Un 
but different special functions. 

Going over from \l/(xi), $(-2i) to the corresponding functions of one variable, \E'(x) 
and $(^), the duality relation ( [4.33|) takes the form 

(1 - zf^^^ (^3^) = -'^ ■ ^(^ ~ ^) *(^) (^■36) 

or, equivalently 

$(x)=C-(-l)^0i^$(i±^), (4.37) 

where we have used that ^{z) transforms to (— z)^$(l — 1/z) under cyclic permutations. 

4.4 Energy spectrum: Exact solution 

4.4.1 Calculation of the energy 

The set of coefficients Uk = Uk{N,q) uniquely defines the eigenfunction ( p.55|) corre- 
sponding to the pair of quantum numbers h = N + ?> and q. Once the eigenfunction is 
known, the corresponding value of the energy £{N, q) can in principle be found by 'brute 
force' as the expectation value of the Hamiltonian. As we show in this section, there 
exists a simpler and much more elegant way to calculate the energy S{N,q) by using 
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the cyclic permutation symmetry. To this end, is proves convenient to work in the dual 
representation. 

The calculation is based on a simple identity Tigs = VHi2'P~^ = VTCl2V^ which 
allows to rewrite the Hamiltonian as 



u^i2 = ^12 + vni^v^ + v^ni^v . 



'•3/2 



(4.38) 



Applying the wave function ^{zi) to the both sides of this relation and using Eqs. ( [4.11| ) 
and ( |3.55| ) we get 



N 



£^{z,) = (1 + eV^ + 9^V) ni^^izi) = Y, «"«n f-^e{n) (l + 9V^ + 9^V) $!i^i 



(12)3, 



Zi) 



n=0 



Here, e{n) denotes the energy of two-particle Hamiltonian ( |3.23|) defined as 

^12 <f = e{n)^fl' , e{n) = 2 [^{n + 2) - ^{2)] . 
Using the explicit expressions ( p.52|) for $)y^ (zj), one can rewrite ( [4.39| ) as 

N 



(4.39) 
(4.40) 



o = J2^-uj;;' 



n=0 



[S-e{n)]^n{z)-e'{~zyip. 



1-i -.(.-i)V^^ 



(4.41) 

which has to be valid for an arbitrary real z. Taking the limit 2; — ^ and using the 
relations (|3.56|) , we get 



£iN,q) 



4(-l)^M, 



N,,-l 




Re 



N 



-^rM„£(n)(2n + 3) 



n=0 



(4.42) 



{N + 2){N + 3) 
Finally, taking into account Eq. ( [4.24|) we obtain the following expression for the energy 

Elo^" «".(?) mn + 2)-^{2)] (2n + 3) 



£{N,q) = 4Re- 



^N 



E„=o^"M?)(2n + 3) 



(4.43) 



To summarize, the recurrence relations ( [4.3|) combined with the expression for the 
energy ( [4.43| ) and the eigenfunctions (|3.55|) provide one with the exact solution to the 
Schrodinger equation for the Hamiltonian 7^3/2- 

An immediate consequence of ( [4.43| ) and the parity property ( [4.8| ) is that 



£iN, q) = £{N, -q), 9iN, q) = l/e{N, -q). 



(4.44) 



Thus, the energy levels corresponding to nonzero values of quantized q are double de- 
generate. 

The resulting spectra of the conserved charge q and the energy £^3/2 are shown in 
Fig. ^ for N < 30. As we are going to argue in Sect. 4.4.3, the eigenvalues form the set 
of trajectories a few of which are shown in Fig. H by solid curves. 
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4.4.2 The exact solution for q = 

The energy and the eigenfunction of the state with g = can be calculated explicitly. 
We recall that this state exists for even N only and its expansion coefficients in terms 
of conformal polynomials are given by (|4.9| ). Their substitution into ( [4.43 ) and ( 4.24| ) 
yields 

S{N,q = 0)=A'${N + 3)+4-fE-Q, e{N,q = 0) = 1. (4.45) 

The curve corresponding to this expression for the energy is shown in Fig. Qd by dots. 
We observe that for even A^ the state with g = is the ground state of the Hamiltonian 
7^3/2. According to (|4.15| ), the corresponding wave function, "^g^Q^Xi), is a completely 



symmetric real function of Xj. Its explicit expression can easily be obtained directly from 



( [4.2| ), without an expansion over the conformal basis. It is straightforward to verify that 
for g = the only solution to ( ^4.2| ) and ( |3.28| ) with the required symmetry is (up to an 
overall normalization factor) 

XiX2X3,^N,q=oiXi,X2,X3) = (4.46) 

= a;i(l - x,)C'Jl,il - 2xi) + 3:2(1 - X2)C'Jl,il - Ix^) + 3:3(1 - X3)&^1^{\ - Ix^), 

where xi + 0:2 + X3 = 1. This translates to 

^^,,=o(x) = 1 - [(1 + x)/2]^+3 - [(1 - x)/2]^+^ (4^47) 

Note that '^N,q=<d{x) does not have zeros on the interval — 1 < a; < 1 and vanishes at the 
end points. 

4.4.3 The Baxter equation, Bethe ansatz and analytic structure of the spec- 
trum 

Numerical solutions shown in Fig. § exhibit remarkable regularity. To understand their 
properties we develop the WKB expansion of the energy £{N,q) and the conserved 
charge q at large N . 

The strategy is in many respects similar to the Bohr's description of the hydrogen 
atom. The Hamiltonian H3/2 describes the system of three particles with the coordinates 
Zi. The scale of the energy is fixed by the conformal spin A^ + 3 which plays the role of 
the inverse Planck constant, h ~ 1/A^, in the corresponding Schrodinger equation. The 
size of quantum fluctuations decreases with A^ and at large A^ the quantum mechanical 
motion of three particles is conflned to their classical trajectories that can be shown 
to have a flnite period. We then quantize the system semiclassically by imposing Bohr- 
Sommerfeld quantization conditions on the periodic classical trajectories. This procedure 
corresponds to the WKB solution of the Schrodinger equation ^{zi) = exp[iN Seik{zi)] 
which, as we will show in the next section, gives a good quantitative description of 
the system. Our aim in this section is develop a physical interpretation of the WKB 
solutions, and to this end we have to introduce some methods of integrable models. 
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The classical analog of the Hainiltonian is obtained by replacing the derivatives by 
the momenta, —idz,, — > Pk, and the commutators by the Poisson brackets. One gets 7^3/2 
as a function of the conserved charges, L^, L^ and Q, each of which describes certain 
modes of the classical motion which will later be quantized giving rise to a complete 
set of quantum numbers. Note that L_ = i{pi + P2 + Ps) is the total momentum of 
the system. The condition ( |3.29| ) then implies that the center-of-mass stays at rest, 
Zi + Z2 + z^ = 0. Similarly, Lq = —i{zipi + Z2P2 + z^ps) + 3 = A^ + 3 generates dilatations 
of the coordinates and its eigenvalue fixes the overall scale of coordinates and momenta. 

The classical motion driven by the conserved charge Q is, however, very nontrivial. 
It generates a collective motion of all the three particles which represents a wave packet 
(or solitonic wave) propagating on the periodic chain with three sites B^ 



z^{t) = Q{Kn + ujt), (4.48) 

where Q{(p) is a 27r periodic function of the argument, and t is the evolution 'time' 
conjugate to the 'Hamiltonian' Q: dzn(t)/dt = {Q,Zn}- uj and k = ^i are the proper 
frequency and the quasimomentum of this wave, respectively, both depending on q and N. 
The periodicity condition Zn+sit) = Znit) leads to quantization of the quasimomentum 
i = integer. The explicit expressions ||32[ can be derived by applying the methods of 
the theory of the finite-gap soliton solutions but are of no relevance for what follows. 
The eigenvalue of the Hamiltonian 7-^3/2 (-^^, <5) defines the energy of the soliton wave 
S = S{N,q). 

Quantization of the charge Q appears as the result of imposing the Bohr-Sommerfeld 



quantization conditions on the periodic classical trajectories ( |4.48| ). To this end, one 
has to identify the corresponding action-angle variables which in turn are constructed 
through the separation of variables. 

The definition of the separated variables for the Hamiltonian 7^3/2 is known thanks 
to the similar construction for the XXX Heisenberg magnet of spin s = — 1 [Q. It 



amounts to the unitary transformation of the operators and the wave function under 
which the original coordinates Zi are replaced by new collective separated coordinates C,i 
and the wave function ^{zi) is transformed into the wave function having a factorized 
dependence on each of new coordinates 

H^^)^Q{^l)Q{^2)^s''-'. (4.49) 

Explicit expressions for the transformation Zi -^ C,i can be found in [^, |28|. The last 
factor in (|4.49| ) carries conformal spin of the state and has a trivial dependence on the 
coordinate ^3. The original 3-body Schrodinger equation for 7^3/2 is translated into the 
Schrodinger equation on the wave function Q(0 and is given by |T6[: 

( AT + 3KiV + 2) ^ ^^j ^^^^ ^ ^^^ ^^^^ ^^^ _^^_ 2^^^^^ ^^_^^^ 

This equation is known as the Baxter equation for the XXX Heisenberg magnet of spin 
s = — 1. In the WKB approach the wave function Q(0 describes the wave function of 
the semiclassically quantized soliton wave in the separated coordinates. 
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We will later show that the Baxter equation is equivalent to the recurrence relations 
([4 .31). Main advantage of considering the Q-function instead of the set of coefficients m„ 
is that it has all intuitive properties of the wave function that one is used to, whereas 
for Un it is difficult to invoke any physical intuition. 

As such, Q(0 should have a finite number of zeros in the classically allowed region 
on the real ^— axis whose position and the total number is determined by the quantum 
numbers N and q. The only 'physical' solution to the Baxter equation ( |4.5CI| ) satisfying 
this condition defines Q(0 as a polynomial in ^ of degree A^ + 3 

N+3 

Q(0 = const X JJ(e-Afc). (4.51) 

fe=i 



Replacing Q(0 in ( 4.5(J| ) by this expression we immediately find that the charge q is 



quantized. Moreover, putting ^ -^ A^ in the Baxter equation it is easy to see that for 
g 7^ one of the roots, Aat+i = AAr+2 = Xn+3 = 0, is three times degenerate and the 
remaining A^ roots satisfy the Bethe equations corresponding to the spin s = — 1 XXX 
Heisenberg magnet 

^^y= n ¥^^' n=l,...,N. (4.52) 

An- V , If, A„- Afc + z 

k=l,k=f=n 



It can be shown [Q that the solutions to the Bethe equations define the set of real roots 
{Afc} which have the properties of the roots of orthogonal polynomials and uniquely 
determine the Q— function ( [4.51J ) as well as the quantized values of the charge 



Q(0 + QH) oT ttA ' 



2Imn(l-Aj^ (4.53) 



lim5^„r'Q(0 ,., 

The explicit expression for the polynomial solution can be found for g = 

■^+i^Ar , Q^^^^ , o^e2 r. ^A^ + 4, -A^ - 1, 1 - 2^ 



Q(0^_^ = Qiv+i(0 = ^"^'(iV + 3)(A^ + 2)e%F2( ' 22 ' ' 1 ) ■ (4.54) 



This expression defines the so-called Hanh orthogonal polynomials [g^ and in the sequel 
we will use some of their properties: 

^"^-^^^ = ~ %n + 3V^ [^"^^(^) + ^"-^^^)] ' 

Q„(±0 = -{TtTin + 2)(n + 1) , (4.55) 

(lnQ„(±z)y = T2t mn + 2) - ^(1)] , 

where the prime denotes a derivative with respect to ^. Since the functions Qn(0 form 
the complete set of orthogonal polynomials we may seek for the general polynomial 
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solution to ( ^4.50| ) for g 7^ as an expansion over the q = solutions 

n=0 ^ n=0 ^ ~'~ ''v "T J 

It is easy to check using (|4.54|) and ( [4.55| ) that thus defined function Q(0 satisfies the 
Baxter equation ( |4.50|) provided that the coefficients m„ satisfy the recurrence relations 
( [4.3D . Thus, the analysis of the recurrence relations (|4.3| ) is equivalent to finding the 
polynomial solutions to the Baxter equation ( [4.50| ). Moreover, comparing the relations 
(4.54) and ( 3.55| ) we observe that the transition to the separated coordinates amounts 



to the replacement i"y9„(2;) -^ Qn(0 i^ the expansion of the wave functions $(-2i) and 
Q(,^), respectively. 

It is worthwhile to note that lengthy expressions for the spectrum of the Hamiltonian 



H3/2, Eqs. ( 4.43 ) and ( 4.24 ), take a remarkably simple form in terms of the Q— function. 
In particular, 

, ^Q'W . Q'(-') „ ^ 2 



To summarize, the Baxter equation ([4.50|) takes the form of a finite-difference Schro- 



dinger equation with the conformal spin A^ + 3 playing the role of the (inverse) Planck 
constant. Applying the standard WKB analysis one can find the asymptotic expressions 
for the solutions corresponding to classical soliton waves propagating on the chain of 
3 particles. The quasimomentum of the soliton is characterized by an integer number 
K = ^i and the proper frequency cj is a (complicated) function of conformal spin A^. 
Changing A^ continuously with i fixed amounts to the adiabatic deformation of the 
soliton solution. This suggests that quantized values of energy and conserved charge q in 
Fig. 2 belong to trajectories parameterized by the integer i defining the quasimomentum 
in ( [4.25| ) and ( [4.48| ). One important property of this deformation which is responsible 



for the analyticity of the trajectories is that it does not destroy the wave packet but 
rather induces the flow of its parameters with A^ known as the Whitham flow ||3^ . The 
precise deflnition of i will be given below. 

4.5 Energy spectrum: WKB expansion 

The WKB solution to the eigenvalue problem for 7^3/2 can be based on the asymptotic 
behavior of the recurrence relations ( [4. 3D at large A^. In this limit it is convenient to 
introduce the scaling variables 

7] = ^Jh{h - 1) = AT + ^ + 0{l/N) , x = r]-\n + ^), q = T^-^q (4.59) 
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such that X takes continuous values on the interval < x < 1. At large N, the 
n— dependent coefficients entering the recurrence relations ([4 .31 ) become functions of 
X, which we define as 



fn ^ V'fix) 



Ur. 



U[X] 



From the definition ( p.33| ) we find the scaling function f{x): 



^(^) = 4^ 



1 - ^ + TT 



x^- 



Aiif 



(4.60) 



(4.61) 



Notice that there is no 0{l/ri) term with our definition of the scaling variables. At 
large N the recurrence relations ( [4.3| ) take the form of the second-order finite difference 
equation 



u{x + f] ^) + u{x — Tj ^) — 2u{x) 



Q 



/(^) 



2 u(x) 



(4.62) 



It has to be supplemented by the boundary conditions 



which can be written as 



u 



1 
2t] 



u 



yrriTv) 



0. 



(4.63) 



The parity property ( |4.8| ) allows to restrict our consideration to positive values of q only. 
It is convenient to interpret Eq. ([4.62|) as a discretized Schrodinger equation with 
?7~^ playing the role of the Planck constant and 2 — q/ f{x) the effective potential. It 
is then clear that the 'wave function' u{x) has different behavior depending on the sign 
of 2 — q/ f{x). The interval of x, on which f{x) < q/2, corresponds to the classically 
forbidden region where \u{x)\ is a monotonous (decreasing or increasing) function of x. 
The crucial observation is that for g > the equation f{x) = q/2 has two real roots, X- 
and x+, on the interval [0, 1]: 

f{x^) = h (4.64) 



for q < l/v27 and none for q > l/v27. In the latter case, |'u(a;)| is a monotonous 
function of x throughout the whole interval < x < 1 and the only way to satisfy the 
boundary conditions ( |4.63| ) is to put u{x) = 0. Therefore, the recurrence relations have 
nontrivial solutions satisfying ( |4.63| ) in the former case only, leading to the constraint on 
possible values of the charge q 



<a< 



1 



27 



27 



(4.65) 



which one readily verifies using Fig. ^ja. For the values of q in this range, u{x) grows 
(decreases) on the interval [0,x_] ([x+,1]) and has a local maximum(s) on the inter- 
val [a;_,x+]. The interval [x_,x+] corresponds to the classically allowed region for the 
Schrodinger equation (|4.62|). 
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4.5.1 Upper part of the spectrum 

We first consider the 'upper' part of the spectrum q -^ l/v27. In this case, x± — >■ 
l/v3 with x+ — x_ = 0{{q — l/v^27)^) and the interval [x_,a;+] shrinks to a point. 
Assuming that u{x) is a smooth function of x on this interval, \u'{x)/u{x)\ <^ f], we 
replace Eq. (|4.62| ) in the leading N ^ oo limit by the second-order differential equation 

1 rf2 



r/^ dz'- 



:U 



(0) 



(z) + 9z'm(°)(z) = --g(^)M(°)(z) , (4.66) 



with z = X — y= and 



u{x) = u^°\z) + 0{r]-^/^) , (4.67) 

g = ^[l + rV'^ + C?(r')] , (4.68) 



which one recognizes as the Schrodinger equation for the harmonic oscillator. Thus, we 
get readily the quantized values of the charge 

g« = -3('£+0, £ = 0,1,... (4.69) 

and the coefficient function 

u^'\z) = H,{^z) exp i-^-^A , (4.70) 

where He{z) are the Hermite polynomials. 

The following comments are in order. The solution (|4.70|) was found under the 



assumption that u{x) is a smooth function, \u'{x)/u{x)\ <^ rj. We verify that it is 
satisfied indeed provided that z <t^ 1 and i <t^ N. For higher excited states, £ ~ A^, we 
are approaching the region g ^ 0, in which the solution is expected to oscillate rapidly, 
cf. ( [4.9| ), and the above approximation does not work. 

The quantized values of the charge, ( [4.69|) and ( |4.68| ), are enumerated by a nonneg- 



ative integer i which counts levels of the harmonic oscillator ( [4.66|) . Using ([4.69|) and 



( [4.68| ) as the definition of the family of curves q = q{N, €) for continuous A^ and discrete £ 
one obtains the trajectories shown in Fig. ^. Namely, the largest values of the quantized 
q for any A^ belong to the same trajectory with £ = 0, the next-to-largest values — to 
the trajectory with £ = 1 and so on. 

The integer d, has a simple interpretation in terms of the solutions ([4.70|) . As ex- 
pected, u{x) oscillates on the interval [x_,x+]. The integer C, counts the number of its 
zeros and the solutions belonging to the same trajectory for different A^ all share this 
number. Recall that considering properties of the exact solutions we have found that 
they are parameterized by a discrete quantum number which is the eigenvalue of the 
cyclic permutation operator ( [4.25[ ). An explicit calculation gives 



e{N, q) = e-^'^^^'"), <j>{N, q) = -{N + i) . (4.71) 
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The approximation ( ^.68| ) can be systematically improved by taking into account 
0{l/ri) corrections to Eq. ( f4.66|) . This allows to evaluate nonleading corrections to the 
spectrum ( [4.68|) and ([4.7CI|) . We obtain 



q{N,i) 



,(1) 



,(2) 



+ ... 



q' 



Q 



(2) 



-3(^+2)' 

2f + 2i 



13 

24' 



,(3) 



,(4) 



9 



(2£^ + 3£^ + 23£+ll), 



1 



31104 



(3840£^ + 7680£^ - 112800^^ - 116640£ - 90899). (4.72) 



Explicit expressions up to 0{ri ^) can be found in [^ ^ . 



4.5.2 Lower part of the spectrum 

In the limit q ^ the classical 'turning points' x± are approaching the end-points 
x„ — i> and x^ ^ 1 where u{x) must vanish. The WKB analysis is not applicable in 
the vicinity of these points and one has to solve the recurrence relation ( [4.3| ) for small n 
and N — n directly, by expanding /„ in powers of n/N and (A^ — n)/N, respectively: 

_^{n + l){n + 2) 



qun = V 



2(2n + 3) 



-[Un+l + U. 



n-1 



U-l 



0. 



qUn = -Tj ^{N - n)[Un+l + Un^l] , Un+1 = 



n<^N, 



N -n<^N. 



(4.73) 



(4.74) 



Quantization of q appears as the condition for these two solutions to match the WKB 
asymptotics in the regions 1 <^n <^ N and 1 -C A^ — n <^ A^, respectively, in which the 
WKB analysis is still applicable. 

For X away from the end-point region, a;_ <^ a; <^ a;_|_ we look for the WKB solutions 
to Eq. ( ^.62| ) in the form 

u{x) = pcos{riS{x)) , (4-75) 

where p is a real normalization factor. Substitution of this ansatz into ( [4.62 ) yields in 
the leading large N limit the following equation on the eikonal phase 

2g 



cos(S"(x)) 



(4.76) 



x{l — x^) 

Solving it at small g ^ 1 we obtain the leading WKB solution for n ^ 1 and N — n^ 1 

as 



Ur.. = pcos ( -n - V? - giVln ^, _ ^, 



pRe 



ju^-t^ 



X 



X" 



-iqN' 



(4.77) 
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with cp being an integration constant. Using Eq. ( |4.24| ) and replacing the sum over n by 
the integral over x one gets 



e-'^ = e, p = {-iy'\V{l-iqN)\-\ (4.78) 

The solution to the recurrence relations for n <^ A^ can be obtained by noticing the 
striking similarity of ( |4.73 ) with the first relation in ( 4.55| ) that describes properties of 



the solution of the Baxter equation for g = 0, so that 

Un = Qn{-qv), (4.79) 

with Q„ as defined in ( |4.54| ). Finally, for (A^ — n) <^ A^ it is easy to verify that the 
solution to (^]7^) is given by 



uj^_^ = z^^-" / dx{l - xY'^'^-'il + x)-*^^-^a;''-". (4.80) 



1 



The three expressions in ( |4.77|) , ([4.79|) and ( [4.801 ) correspond to the solution of the the 



recurrence relation in the three different regions which overlap however, for 1 <^ n <^ A^ 
and 1 <t: N - n <^ N. Requiring that ( [4.77|) can be sewed with ( [4.79|) for 1 < n < A^ 



and with ( [4.80| ) for 1 <^ A^ — n -C A^, we find the quantization condition on q 



qN~'\nN-argT(l + iqN~') + 0(l/N) = -{N-2i), (4.81) 

6 

where i is an integer. This result is valid to 0{1/N) accuracy for small q/N^ <^ 1 and 



can significantly be improved by taking into account nonleading corrections to ( |4.81| ). 
In this way one gets 

7if{q, N)/3 = q*r]lnr]- argr(l + ir]q*) + |- - —riq*' + Oim*') = T^iN - 2£), (4.82) 

or] b o 

where 

q* = q{l + 2f). 

For given i, the quantized values of q belong to the i—th trajectory q = q{N, i) which 
depends analytically on A^. It follows from ( [4.82| ) that the function q{N, i) has the 
reflection symmetry 

q{N,i) = -q{N,N-i), (4.83) 

which maps positive values of q on the i—th trajectory into the negative q on the (A^ — 
£)— trajectory. The i—th trajectory crosses the zero g = at even N = 2i and rises 
towards larger values of q corresponding to the 'upper' part of the spectrum. 

To illustrate this property, we evaluate the function f{q, N) for the exact numerical 
values of q belonging to the same trajectory shown in Fig. ^ja, and plot it against A^ as 
shown in Fig. |^. It is seen that the linear behavior in A^ continues to the upper part of 
the spectrum where Eq. ( |4.82| ) can be matched with the WKB expansion in ( |4.72| ). In 
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Figure 3: The trajectories for q as given by Eq. ( |4.82| ). The crosses give the values of the 
function f{q,N) calculated using the exact eigenvalues q. 



this way we can check that the definitions of i in Eqs. ( 4.72|) and (^4.82 ) do match each 
other and describe the same trajectory. 

Solving Eq. ( [4. 821 ) for small N — 2i and large r] one gets 



q/r]'^ 



IT 



-(iV-2£) + C((£/lnr/)3) 



(4.84) 



so that a few lowest eigenvalues of Q are of order O^r]"^/ Inrj). 

4.5.3 Asymptotic expansion of the energy 

Let us use the WKB solutions to the recurrence relations to obtain the asymptotic 
expressions for the energy. 

For the 'lower' part of the spectrum we substitute ( [4.77] ) into the exact expression 
for the energy, Eq. ( [4.43| ), to get after the integration 

£^ = A ln{N + 3) - 6 + 67ij + 2Re ^^(1 + iN'^q) . (4.85) 

This expression is valid for q = 0{N'^) up to corrections suppressed by powers of 1/A^ 
and, in particular, for g = it reproduces the exact result ( [4.45|) . 

A more accurate and general expression can be obtained [^ Q by asymptotic 
expansion of the Baxter equation: 



£g = 2 In 2 - 6 + 67i5 + 2Re ^ V(l + ^rf^k) + ^(r/-^). 



(4.86) 



fc=i 
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where 5k are defined as roots of tlie following cubic equation: 

(5fc - g = 



n 



(4.87) 



and q satisfies the condition ( [4.65| ). It is easy to see that for q belonging to the interval 
( [4.65| ) all roots 5k are real. The expression in ( [4.86| ) is valid with high accuracy for the 
whole spectrum. For q = 0{N'^) both expressions, ( [4.85| ) and (|4.86| ), coincide. 
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Figure 4: The dependence of the energy £ on the charge q for A^ = 30. The soUd curve is 
calculated using Eq. (4.86) and the exact values of energy for quantized q are shown by crosses. 



The resulting dependence of the energy Sq on the charge q for A^ = 30 is shown in 
Fig. §. We find from ( 4.85 ) and ( [4.86 ) that the energy is quadratic in q close to the 
origin g = 

2C(3) o 



£g = A\n{N + 3) - Q + A-fE + 



N^ 



-q 



(4.88) 



with C(3) = 1.20205690, whereas at large q = 0{N^) the asymptotic behavior of the 
energy is given by 

£g = 2\nq-6 + 6-fE + 0{N-^). (4.89) 

We would like to stress that the expression ( [4.86|) defines the dependence of the Hamil- 
tonian on the conserved charges 7^3/2 = 'Hz/2{Q, L"^) for large eigenvalues of L^. To find 
the spectrum of the energy £ one should replace q in (|4.86| ) by their quantized values. 

Calculating the quantized values of the energy £q we find that each trajectory g(£, A^) 
is mapped into the corresponding trajectory for the energy £{l, N) as shown in Fig. 0b. 
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In particular, the i—th trajectory starts a.t N = i, approaches the 'Fermi surface' Sq=o 
a.t N = 2i, gets repelled from it and monotonously grows to infinity at large N. The 
corresponding asymptotic expression for the energy reads p9|, |30| : 

S(N,i) = 6\nr] - 3\n3 - 6 + 6-fE --i'2i+ I) -i, (5 f + 5i -7/6) 

r] 1]^ 

- -J— (464^2 + 696f - 802£ - 517) + . . . (4.90) 

at large N ^ d, and 

£:(iV,^) = 41n(iV + 3)+47i;-6+ J"y^ (iV-2/)^ (4.91) 

18 m [rje^'^) 

in the vicinity of A^ = 2L To find the behavior around N = i one has to use Eqs. ( |4.83| ) 
and (|4.44| ) to get 

S{N,i) = S{N,N-£) (4.92) 

and substitue the expression in the r.h.s. by ( [4.90| ) with i replaced hj N — i. 

The relations ( |4.90|) , ( |4.91|) and ( [4.92|) define the asymptotic expansion for the energy 
levels of the Hamiltonian 7^3/2 parameterized by the integer i. We observe that for given 
A^ the distribution of levels is different in the lower, i ~ N/2 or equivalently q -^ 0, and 
the upper part of the spectrum, i <ti N oi q ^ l/v^27- Using ( [4.90| ) and ( [4.91| ) we find 
the corresponding level spacings as 

6£{N, i) '=' O (-\-\ , 5£{N, i) '^'i^ O ("-] . (4.93) 

4.6 Analytical continuation and the parton model 

Each polynomial eigenstate of the Hamiltonian 7^3/2 corresponds to a multiplicatively 
renormalizable local operator and thus an independent nonperturbative parameter in the 
distribution amplitude ( |2.19| ). If, as usually assumed, the sum in ( |2.19| ) is uniformly con- 



vergent pointwise in Xj, then the baryon distribution amplitude is restored uniquely from 
this expansion. The assumption of uniform convergence ensures that the distribution 
function vanishes as xxx^X'i at the end points x, — > and implies that the nonpertur- 
bative reduced matrix elements decrease sufficiently fast for large conformal spins. If 
the initial condition to the Brodsky-Lepage evolution equation (distribution amplitude 
at a low scale) decreases for Xj — ;> at a slower rate, then the series in ( |2.19| ) diverges 
close to the end points and the scaling behavior has to be defined by a (infinite) resum- 
mation of the dominant contributions of large conformal spins. This resummation can 
be performed in the standard way by replacing an infinite sum over A^ by an integral 
over complex A^. To this end the analytic continuation in A^ becomes necessary and, in 
particular, the anomalous dimensions ^fq ought to be analytical functions of A^. It is this 
situation that occurs in the study of the x -^ 1 and x — > limits of an inclusive process 
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with exchange of baryon quantum numbers in which case forward matrix elements of 
baryon operators contribute and the expansion in moments leads to the expansion of the 
corresponding generalized parton distribution in derivatives of the 5-function at x = 0. 
As familiar from studies of deep inelastic scattering, restoration of parton distribu- 
tions from known values of the moments calculated within the framework of the operator 
product expansion involves, first, an analytic continuation of the anomalous dimensions 
from integer positive A^ corresponding to spins of composite operators into the complex 
A^ plane and, second, decomposition of the distribution amplitude 0(xj; /x) for arbitrary 
Xi into irreducible components having an appropriate analytical (spectral) properties 
and admiting the parton model interpretation. This procedure is well understood for 
leading twist parton distributions, see e.g. P3[, but, to our knowledge, has never been 



discussed for three-particle distributions. 

Mathematically, the first task consists of defining anomalous dimensions as analytic 
functions of A^ such that their values at positive integer A^ are given by the eigenvalues of 
the evolution kernel 7^3/2 and the asymptotics at infinity is such that 7Ar ~ exp(— (5|A"|) 
with 5 < n ^E\. The anomalous dimensions of the different components of the dis- 



tribution amplitude do not necessary coincide and it is known from the studies of the 
deep inelastic scattering that one may need to consider analytic continuation from odd 
and even A^ separately, which in general correspond to contributions of operators with 
different parity. 

The set of trajectories shown in Fig. 2 presents a legitimate analytic continuation 
and reflects the highly nontrivial analytic structure of the integrable model. This set 
is complimentary to a simpler and more general analytic continuation corresponding to 
ordering of the anomalous dimensions from below, see Fig. ^. The trajectories shown in 
Fig. I^a are copied from Fig. ^. They are enumerated by an integer number i defined 
in ( [4.69|) the physical interpretation of which is discussed above at length. Going over 



to the trajectories shown in Fig. ^d corresponds to the rearrangement of eigenvalues 
according to a different integer number i which is related to £ by a formal substitution 

e = [N/2] - i, (4.94) 

where [A^/2] denotes the integer part of A^/2. The expression on the right hand side 
of Eqs. ( |OTl ), (|]82|) then becomes 7r/3(£ + (5iv/2) where 6n = and 6n = 1 for even 
and odd values of A^, respectively. The assignment of the eigenvalues to trajectories 
becomes, therefore, different for odd and even A^. The three trajectories shown in Fig. ^3 
correspond to £ = 0, 2, 7 and correspond to the analytic continuation from even A^. In 
particular, the £ = trajectory going through the lowest eigenvalues at even A^ is given 
by Eq. ( [4.45| ). Note that the corresponding eigenstates all have positive parity. The two 
(degenerate) lowest trajectories for odd A^ formally correspond to i = and £ = —1 
and can further be rearranged in contributions of definite parity. Note that, in contrast 
to Fig. ^, each trajectory in Fig. ^d corresponds to a fixed eigenvalue 6 of the cyclic 
permutation operator. Another important difference is that each trajectory in Fig. ^ 
behaves as ~ 61nA^ at A^ ^ oo while each trajectory in Fig. ^ grows as ~ 41nA^ 
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Figure 5: Two different analytic continuations of eigenvalues for the integrable Haniiltonian 
7^3/2, see text. 

only. The upper boundary ~ 6 In A^ (shown by dots) arises in this case because new 
trajectories are being built on the top of the spectrum starting at each integer A^. In 
both cases the asymptotics of the anomalous dimension at large N does not exhibit an 
exponential growth and an analytical continuation to complex N is unique. 

We emphasize that both sets of trajectories define legitimate analytic continuations 
and the one shown in Fig. ^ is made possible by an additional 'hidden' symmetry on 
the integrable Hamiltonian. The choice between them is defined by the process in which 
the baryon distributions are measured or, equivalently, by the way in which the end- 
point region in (xi, X2, X3)— space is approached. It remains to be studied which analytic 
continuation ensures the true asymptotic behaviour of the distribution amplitudes at the 
end points. This question goes beyond the tasks of the present paper. 



4.7 Eigenfunctions 

In this section we find an explicit expression for the eigenfunctions \E'(x) in the large N 
limit. Requiring Q\l/ = g\l/ yields a third-order differential equation 



(h - l)(/i - 2)^'(x) - 2{h - 2)x^"(x) - (1 - x2)^'"(a;) = -2iq ■ 



-^(x) 



(4.95) 



l-a;2 

where h = A^ + 3, which is symmetric under the transformations corresponding to cyclic 
permutations Xi ^ X2 ^ X3 -^ Xi 

h 



X 



V 



a; — 3 

x + r 



^(x) 



V 



-v. 



x + l 



^ 



X — 3 

x + l 



(4.96) 



45 



V 



V 



Note that the cychc permutation symmetry maps the interval [—1, 1] -^ [— c>o, — IJ 
[1,00]. It is sufficient, therefore, to consider the region — 1 < x < 1 only since the 
function \E' outside this interval can be recovered by the transformation ( |4.96| ). 



4.7.1 WKB solution 



Eq. [4.95| can be solved at large values of r^ = \/h{h — 1) by the WKB expansion. To this 
end we write 



^(x) 



exp 



r7 5o(a;) + 5i(x) + 0(l/r7) 



(4.97) 



with 5*0, 5*1, . . . being 77- independent. We are going to use that values of the quantized 
charge q = q/iff are smaller than (f < 1/27, see ( |4.65 ). This allows to expand the 
functions 5*0, 5*1, . . . in powers of g <^ 1. 

In particular, substituting the WKB ansatz ( [4.97| ) into the differential equation ( |4.95| ), 
one gets in the leading 77 ^ 00 limit the following equation for S'q = dxSo{x): 



(1 - x')S',ix)[l - (1 + x)S',] [1 + (1 - x)S',] = -2i q. 



(4.98) 



This cubic algebraic equation has three independent solutions related to each other by 
the symmetry transformation ([4.96|) . Therefore, it is sufficient to consider only one of 
them, the one which vanishes as g ^ 0. The first few terms of its expansion in powers 
of q are given by 



S',{x) 



2iq 



^q^ X 



1 



x^ 



{l-x'^f 



+ 8iq' 



'^^^'^ +I60q'x^±^,+O{f). (4.99) 



1 — X 



2^3 



1-X 



2^4 



Integrating this result and substituting 6*0(0;) into ( [4.97| ) one gets the leading WKB 
asymptotics 



"^(x) 



1 + x 



1 



X 



-iriq 



exp 



4:1] q 



l-a;2 



(4.100) 



where the 0{q^) and 0{q'^) terms are omitted for simplicity. To find the first nonleading 
correction to this expression, one further expands the differential equation and keeps the 
0{l/ri) terms. In this way, one finds the expression for S'i{x) in terms of Sq{x) and its 
higher derivatives that one integrates to get 



Si{x) 



const 



iiqx iq ^ 1 + x 
o + ^ In 



1 - x^ ■ 2 """ 1 - X 
Finally, one obtains the following WKB approximation for the eigenfunction 



(4.101) 



^ 



WKB, 



X 



const X 



1 



X 



X 



-irjq 



exp 



4:T]q 



^2 



Aiqx 



x^ 



x^ 



+ o{m'') + o{i/ti) 



.(4.102) 



The constant fixes the normalization of \^ and is otherwise arbitrary. Applying the 
symmetry transformation ( [4.96| ) to this expression one can construct the two remaining 
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Figure 6: Exact (solid lines) and WKB (dashed lines) eigenfunctions 'I'(x) for A^ = 60. The 
figures in the first and the second row correspond to the maximum value of q (alias energy) 
and minimum nonzero value of q, respectively. The left and the right figures show real and 
imaginary parts of ^(x), respectively.. 

fundamental solutions to the differential equation. The general solution is given by 
their linear combination with two arbitrary constants. The latter can be fixed from the 
requirement for \E' to be an eigenfunction of the permutation operator. In this way, one 
finds the final expression for the eigenfunction as 



'^(x) 



^ 



WKB 



+ ^(-1) 



1 — X 



x + 1 



^ 



WKB 



WKB 



^ 



X + 3 
1 — X 



X 



1 + x 



(4.103) 



This approximation is valid for large i] and for all real x except in the vicinity of the 
singular points x = — l,l,oo. Note that the two added terms in (|4.103|) represent 
'quantum' corrections to the WKB solution in ( |4.102| ) in the region —1 < x < 1, which 
are not seen to all orders in the l/rj expansion. 



4.7.2 Resummation of leading corrections 

As seen from the above, the expansion in powers of q actually proves to be the expansion 
in 4g/(l — x^) and is compromised close to the end points. One can improve the WKB 
expansion in the previous subsection by making a resummation of the leading singular 
[g/(l — xy]'^ terms to all orders. To this end, we consider the limit a; — >■ 1 of the 
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differential equation Eq. ( |4.95| ): 



S',[l-2S'^] = -iq/{l-x), 



where from 



On 



1-W1 + 



■iiq 



X 



(4.104) 
(4.105) 



Integrating this relation and adding a similar contribution from x ^ — 1 we get 



sr{x) 



1 + x 



-iqhi 



1 



1 



X 



2zgln 



1 



liq/{l 



X 



1 + v/l-82g/(l + a;) 

1 + ^1 + 8ig/(l - x) 

1 — X 



1- ^/l + 8iq/{l-x)] . (4.106) 



4 V ^ ^' ^ V 4 

Adding the less singular terms and collecting everything, we get the leading-order re- 
summed WKB eigenfunction: 



* 



WKB, 



X] 



const X 



X exp 



1 



X 



-im{i+2f) 



1 — X 
1 + x 



l + ^/l-8iq/{l + x) 
1 + V^l + 8zg/(l - x) 



-2irjq 



-7] 



X exp < —r] 



4 
1 — X 



^/l-8iq/il + x) 



4 



X exp 



4iq{l + riq^)x 
l-x2 



- 8iq/{\ — X 
40r/g^ 



'1 — X 



2^2 



+ 0[r]q\l - x' 



(4.107) 



which has to be inserted in (|4.103|) and presents our final result. 

The numerical comparison of the exact and the WKB eigenf unctions in presented in 
Fig. ^ for A^ = 60. This large value of N is chosen to illustrate that the eigenfunctions 
corresponding to large energy eigenvalues (the two upper figures) have a typical wave 
packet structure: The size of the packet is of order 1/yN and the oscillation frequency 
of order ~ 1/A^. The eigenfunctions corresponding to lowest eigenvalues are, on the 
contrary, smooth functions (the two lower figures) for which the WKB approximation 
works very well. 



5 Helicity A = 1/2 distribution amplitudes 

The scale dependence of the A = 1/2 distribution amplitudes is driven by the Hamiltonian 
7Yi/2 defined in ( |2.26|) , which differs from 7^3/2 by the two terms corresponding to gluon 
exchange between quarks of opposite chirality, see Fig. |I|: 

^1/2=^3/2 + 1^, V ' ^ 



-^12 



7'2 

-^23 



(5.1) 
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Figure 7: Matrix elements of the exchange interaction \{'^qfii\\l/L\2\^qU))\ evaluated between 
the 7^3/2 eigenstates labeled by the integer (.,(.' = Q, .. .N + 1 which enumerates the quantized 
values of q and q' from the above, see Sect. 4.5, for N = 30. The picture to the right shows 
the contour plot of the 3-dimensional plot. 



The spectrum of eigenvalues of 7ii/2 corresponds to the spectrum of anomalous dimen- 
sions in the evolution equation. The Hamiltonian 7ii/2 is not integrable and the corre- 
sponding eigenproblem cannot be solved exactly. For a given A^, the spectrum and the 
eigenf unctions can most efficiently be calculated by a numerical diagonalization of the 
mixing matrix for the additional — exchange interaction — terms evaluated in the basis 
of the exact eigenstates of 7^3/2- The reason is that this matrix is strongly peaked at the 
diagonal: Matrix elements {'^ q'\\ / L1^^^'^ q) between the 7^3/2 eigenstates labeled by the 
values of the conserved charge q and q' decrease rapidly with \q — q'\, see Fig. ^ In con- 
trast to 7^3/2, the Hamiltonian 7ii/2 is not invariant under cyclic permutations but, still, 
is symmetric under the interchange of the first and the third quarks, [7^1/2,7^13] = 0. 
This allows to choose its eigenstates to have definite parity with respect to the Pis 
permutations 



7^i/2^S(x.)=^;;2^(iV,£)vl/W(x. 



Vu^^tN^^i 



±^S(^. 



(5.2) 



Here, A^ refers to the total number of derivatives and d, numerates the energy eigenstates. 
It is therefore natural to decompose the eigenfunctions '^'^^\xi) over the basis of eigen- 



states of 7-^3/2 with definite parityQ \l/^^(a;j) and g > 0, defined as in ( |4.2(J| ). Using the 



■ q,N\^' 
-2-02 



identity L^2 + -^23 — "^ -^13 ^ + "^-^13 "^ ^^"^ taking into account that the states \l/g 

^•^Throughout this section we define parity with respect to permutation of the first and the third 
quark, instead of the first and the second quark in Sect. 4. To account for this change, we also use the 



basis of the 'permuted' 7^3/2 eigenstates 7^23^^, with ^it- defined in (4.20). We usually omit 7^23 and 



use the same notation ^* to simplify the presentation. 
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Figure 8: The spectrum of eigenvalues for the Hamiltonian TCi2- The lines of the largest and 
the smallest eigenvalues of 7^32 are indicated by dots for comparison. 



diagonalize the cyclic permutation operator V we obtain 



,(±) 



<t^i^iO = -2/Ap^E 



m=0 



^ fm^Um{q')Um{q) 

(m+ l)(m + 2) 



cos 



-0,o±(-ircos(0,+</>,o], 

(5.3) 
where M„(g) correspond to the expansion coefficients of the 7^3/2 eigenfunctions in the 



m 



(31)2 

n,N 



basis and 



'M^\m<^^^l) 



The factor 



■q',N\ 



N 



■q,Nl 



0. 



^/ = 5Z fm^Um{q)Um{q) 



(5.4) 



(5.5) 



m=0 



comes from the normahzation condition for the states, which we assume in this section 



tobe(vl/W|vl/W) 



1. 



The exphcit calculation gives the spectrum shown in Fig. ^. The lines of the largest 
and the smallest eigenvalues of 7^3/2 are indicated by dots for comparison. 

As seen from the figure, the spectra of ?^i/2 and 7^3/2 are very similar in the up- 
per part, for larger eigenvalues, and at the same time the two lowest levels of the 
7Yi/2 Hamiltonian appear to be special and 'dive' considerably below the line of low- 
est eigenvalues of 7-^3/2, given by Eq. ( [4.45| ). Our goal in this section is to explain 
this structure and to get the quantitative description of the 7^i/2 spectrum in the 
large N limit. Note that at A^ ^ 00 the spectrum of 7^3/2 becomes very dense, see 
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Figure 9: The flow of energy eigenvalues for the Hamiltonian 'H(e) for A^ = 29 and A^ = 30, 
see text. The solid and the dash-dotted curves show the parity-even and parity-odd levels, 
respectively. The two vertical dashed lines indicate 7^3/2 = Ti-ie = 0) and TC1/2 = 'H{e = 1), 
respectively. The horizontal dotted line shows position of the 'ground state' given by Eq. (|4.45| ). 



Eqs. ( [4.93|) , and approaches a continuos spectrum inside the band of the width ~ 2 IniV: 
4 In A^ — 6 < £^i2 < 6 In A^ — 6 — 3 In 3. We will demonstrate that exactly the same band 
of the continuos spectrum is formed for the 7^1/2- Inside the band, the distribution of 
levels is perturbed by corrections at most 0{1/N'^) and 0{l/\r? N) at the upper and 
the lower boundary, respectively. In addition, the two lowest eigenstates of 7ii/2 (one 
for each parity) fall below the 'Fermi surface' and are separated from the bottom of 
the band by a finite constant. Existence of such a 'mass gap' presents our main result 
in this section and its formation will be interpreted as due to binding of quarks with 
opposite chirality by the exchange interaction and formation of scalar diquarks. The 
eigenfunctions of the 'bound states' and the value of the 'mass gap' will be estimated. 

To visualize both the similarities and the differences between the spectra of 7^3/2 and 
?ii/2 and to trace formation of the 'mass gap' for 7^i/2, it proves convenient to introduce 
a somewhat more general Hamiltonian 



'H{t) = n,/2 + tV. 



(5.6) 



with e being a new coupling constant. 7Y(e = 0) reproduces 7^3/2 whereas 7-^(e = 1) 
coincides with the Hamiltonian 7ii/2. Thus, the spectra of 7^3/2 and 7ii/2 are related to 
each other through the flow of the energy levels of 7i(e) from e = to e = 1, see Fig. ^ 
Note that for e 7^ the Hamiltonian 7i(e) is neither integrable, nor cyclic symmetric. 
It is still invariant under conformal transformations and under permutations of the flrst 
and the third quarks, ['H(e), P13] = 0, but the degeneracy between parity-odd and parity- 
even eigenstates is lifted and, in fact, the flow of levels with different parity is completely 
independent from one another. 

The spectra in Fig. exhibit the following characteristic features: 
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In the upper part of the spectrum the effect of e— proportional terms on the spec- 
trum of the 'unperturbed' Hamiltonian ?i(e = 0) = 7^3/2 is very mild. While at 
e = the energy levels are double degenerate, their splitting at e 7^ remains 
(exponentially, as we will argue) small for large N. 

For e > 0, the two lowest levels are decoupled from the rest of the spectrum and 
fall off with e almost linearly. For e < 0, the levels with different parity start to 
cross each other, whereas the flow of the levels with the same parity follows the 
pattern well known in quantum mechanics as the 'repulsion of levels' 
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This structure suggests that the difference eV = Hi^e) — H{e = 0) can be considered 
as a perturbation for most of the levels, but not for the few lowest ones (for large 
A^). To formalize the argument, one has to evaluate the matrix elements in ( p.3|) and 
compare them with the energy splittings for the 'unperturbed' Hamiltonian. The explicit 
calculation (see below) gives 



'? 



(n^^l^l^.iv 



(±)\ g.g^o 



o 



InN 



'^f^ 

"?': 






(±), q,q'^_l./V27 



o 



1 

iV2 



(5.7) 



Comparing this result with the level spacings in Eq. ( [4.93| ) , we conclude that the pertur- 
bation theory in V is justified for large N (or for small e <^ 1 and arbitrary A^) for the 
upper part of the spectrum, while several (of order ~ In A^ as we will find) lowest energy 
eigenstates are affected strongly and have to be rediagonalized (unless e <^ 1/lniV). In 
the sequel, we are going to consider the two different regions separately in more detail. 



5.1 Upper part of the spectrum 

Eq. (|5.3|) becomes, for diagonal transitions q' = q 



(^)1S = (^S^^|V^I^?^) 



N 



-2A/;E 



/, 



m m 



(9) 



m=0 



'm+ l)(m + 2) 



[l±(-l)'-cos(20,)] . (5.8) 



Each term in the sum is explicitly positive so that the matrix element {V) j(f is always 
negative meaning that for all levels £3/2 — ^1/2 > 0. According to ( [4.70| ), the coefficients 
Um{q) are smooth functions of the scaling variable x = m/N, peaked around x = l/y/S 
and rapidly decreasing outside the region {x — l/-\/3)^ < r]~^/3. Splitting the sum in 
( ^.8|) into contributions of even and odd m one finds that the terms proportional to 
cos{2(f)q) tend to cancel each other and their total contribution is approximately given 
by the sum of two boundary terms. This contribution is negligible compared to the sum 
of the phase-independent terms, but at the same time it defines the splitting between 
energy levels with the different parity 



ci+) 






i-irf-'ul{q) 



(m+ l)(m + 2) 



(5.9) 
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Here m^i^/N <^ 1 and 1 — nimax/N ^ 1 define the interval, mmin < m < mmax, on 
which the WKB expansion ( [4.70| ) is applicable. Applying similar arguments to the sum 
entering the matrix element 



N 



= (^.1^-,) - Y.^-irf^'^'^rniQ) , (5.10) 

we conclude that the values of (— l)™/^^^ have to coincide at the end points so that 

^'/^ ^ ^"^ ~ (m.„.,+ l)(m„„., + 2) ~ " ^^P'-"' ■ f^'"' 

Thus, the splitting between the energy levels of different parity is governed by the tail 
of the wave function ( [4.7U| ) and, as a consequence, is exponentially small at large N. 

Neglecting exponentially small terms and, in particular, the level splitting £y^ ^^1/2 5 
we can replace the sum in (|5.8| ) by the integral over the scaling variable x = m/N and 
substitute the coefficients Um by their WKB expansion ( [4. 671) with the leading term given 



by ( |4.70|) . Substituting this expansion into (|5.5| ) and ( |5.8|) one can calculate the leading 



and the next-to-leading corrections to the energy ( p.8| ) a^ 



s[f^{i)-£s/2ii)^{V),,N = -^ 



\ I 2(2f + l) 



+ 0{W), (5.12) 



r] 
verifying the estimate in (|5.7| ). We remind that r] = \/{N + S){N + 2). 

5.2 Lower part of the spectrum 

The analysis of the low part of the spectrum is considerably more involved. 

We start with calculation of the normalization factor Afq defined in ( ^.5|) . Assuming 
q = q/r]^ <^ 1 one can use the WKB approximation ( [4.77| ) for Um{q) to get 

Ar,Hr(l-.?/V)|^('|'l//J .|,M„-(,e™)(-||^)\ (5.13) 



n=0 



This expression defines Mq to be an exponentially decreasing function of \q\N . Because 
of this factor, most of the elements of the (A^ + 1) x (A^ + 1) matrix {^^q^\v\^^q^^) 
(jTS) are very small, see Fig. |^. It is clear that the off-diagonal matrix elements of this 



matrix are those responsible for the mixing of different energy levels of the 'unperturbed' 
Hamiltonian and, therefore, strong mixing can only occur if qN, q'N ~ 0{1). Using the 
large- A^ approximate expressions for quantized q ( [4.84|) we may expect that the number 
of such levels /cmax is of order In 77. 

^''To calculate the next-to- leading correction to this expression one needs to know the O (77^^/^)— term 
in the WKB expansion ( 4.67 ) oiu{x). It can be obtained by substituting ( 4.67 ) into ( 4.62| ) and comparing 



the coefficients in front of the nonleading power of I/77. 
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,(±) 



(±)\ 



To calculate the matrix elements (\E'g |V|\E'g ) given by ( |5.3| ) we observe that, unlike 
in the calculation of A/'g, the sum over m is saturated by the contribution of the first few 
terms, m <^ N . Hence Um{q) can be replaced at small qN by U2m = (—1)™ + 0{qN)^ 
and U2m+i = 0{qN), corresponding to the solutions to the recurrence relations ([4.73| ). 
This gives 



{^f\V\^?^) = -2U, Yl 



J n 



n=even 



(n + l)(n + 2) 



cos 



(0g - (t>q')± COS(0g + 0g/)] 



TT 



— — In ^{rie^'^) [cos{(j)g - (j)qf) ± cos{(f)q + 0^/)] 



(5.14) 



The g— dependence enters this expression through the phases (f)q and 0' which take 
quantized values defined in ( [4.12| ). It is easy to see that the possible values of cos{(f)q± 



) 
are 1 and —1/2 depending on whether the phases (pg and 0^/ coincide. For the present 

purpose it turns out to be more convenient to use the basis of eigenstates "^g with fixed 

q rather than fixed parity "^g ~ \E'q ± \E'_q and write (^.14 ) in matrix form, introducing 

an integer k = [N/2] — i, k = 0, ±1, ±2, ... to numerate quantized values of q starting 

from the ones with the lowest absolute value. We get: 



/ 1 



{^g,\eV\^g) = -gA 



k'k 



^k'k 



2 2 2 

I -I -I 1 



\ 



(5.15) 



where, to our accuracy. 



q 



k- 



Tirj 



and the dependence on e and rj is absorbed in the 'effective coupling' 



(5.16) 



g 



en 



91n(?7eT'E) 



(5.17) 



The approximation in ( |5.15| ) is justified for \k\ -C /cmax = 0{lnr]). 



Imagine, for a moment, that all the participating levels, \k\ < fcmax of the 'unper- 
turbed' Hamiltonian 7^3/2 were degenerate, i.e. their energy splitting negligible compared 
to the interaction in ( ^.15 ). The true energy eigenstates would coincide then with the 
eigenstates of the mixing matrix A^^/ (of the size fcmax), and the corresponding eigenval- 
ues would define the energies. Remarkably, the spectrum of Akk' is extremely simple. 
One can easily convince oneself that Afc^/ has two and only two nonzero eigenvalues which 
both are equal to fcmax- Remembering that g ^ 1/ Inr] and fcmax ~ Inrj, this implies that 
two energy levels will get shifted by the finite amount gk^^;^ = 0(\n^r]) while all the 
other ones remain exactly degenerate to this accuracy. Since (|5.15| ) was derived up to 
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corrections of order ~ l/ln^ r], this implies that the true energy shift for all levels apart 
from the lowest two ones is at most ~ 1/ In^ rj. This simple heuristic observation explains 
the pattern observed in Figs. H, p. 

The real situation is certainly much more complicated. The splitting between lowest 
energy levels of 7^3/2 cannot be neglected, and, in fact, it is precisely this energy splitting 
which determines the number of lowest states /Cmax that can effectively be considered as 
degeneratep^ and the precise value of the 'mass gap'. For small values of qN we can 
calculate the energies S3/2 using the asymptotic expression in ([4.91|) that we rewrite as 



^3/2 (g) = ^0 + 



2C(3)vr^ 



k + 



6 



N 



+ 0{{qNY 



qN 



TX 



6 \\i{rje 



he' 



(5.18) 



where k is defined as above, Eq is the ground state energy given by ( [4.45| ) and ^at = , 1 
for even and odd A^, respectively. Combining together ( p.l8| ) and ( |5.15| ) we obtain that 
in the lowest part of the spectrum, corresponding to qN , q'N < 1, the Hamiltonian 7i(e) 
can be represented by the following matrix: 



2m 



3 V 2 



Sk'k — At' 



k'k 



(5.19) 



IT 



m 



9C(3) 



eln(?7e'^^) 



The corresponding eigenvalue problem is solved in Appendix B. The idea of the 
solution is to interpret the integer /c as a discrete momentum variable. Then, the ex- 
pression in ( p.l9| ) can be considered as an effective Hamiltonian for the low 'frequency' 
1^1 < ^max = 0{\n.rj) modes of 7^(e) and the two terms in the r.h.s. of ( |5.19| ) can be 
identified as the kinetic energy and the periodic potential for a particle on a line. The 
corresponding wave functions in configuration space correspond to Bloch-Floquet waves 
and the resulting Schrodinger equation turns out to be a generalization of the famous 
Kroning-Penney model of a single particle in a periodic (5-function potential^. The 
solution then follows the classical procedure PB| . 

One has to keep in mind, however, that the effective Hamiltonian in ( |5.19| ) presents an 
approximation to 'H{t) up to corrections of order 0{qN) and one has to check whether 
values of qN are small on the solutions. It is possible to show that this condition 
is indeed satisfied for small e <^ 1 (and this is the place where introducing e as a new 
parameter starts to play a role), see Appendix B. From this analysis we obtain, therefore, 
a quantitative description of the spectrum of ?i(e) in the region 



l/ln?7 < e< 1 



(5.20) 



^^So that for the state with the number kmax this energy splitting is of order of the largest eigenvalue 
of the matrix A. 

^^Inclusion of a few first corrections in qN to the matrix A effectively amounts to smearing of the 
(5-function potentials. 



55 



£k = Eo+ ::' _, , (5.21) 



where the lower bound comes from the condition that the interaction eV is sufficiently 
strong to excite many levels. 

In agreement with the heuristic argument given above, we find two bound states and 
the continuum spectrum. The levels in the continuum are e-independent, in the small-e 
limit, and are given by 

81n^(r/e'Ts) 

with 1 < A; < In 77. The two bound states are degenerate up to l/?7^— corrections that 
we neglected from the beginning and the binding energy (which we identify as a 'mass 
gap') is given by 

72C(3) 
Comparing this expression with the nonzero eigenvalue of the perturbaton in ( |5.15| ) 



A(e) = Abound - ^0 = -e'^r^. (5.22) 



5'^max, we calculate the number of excited levels of the unperturbed Hamiltonian 7^3/2 
as fcmax = CTT^ In 77/(8^(3)), which is in agreement with our expectations. 

The low-frequency part of the wave functions of the two bound states coincide for 
ln?7 -^ 00 with the two lowest energy eigenfunctions of V, that are given by the parity 
even and parity odd combinations of the basis functions ( |3.38D diagonalizing the Casimir 
operators L^2 and L"^^: 



^£L(^o 



¥^lUx,) ± vi/gLo(^0 



+ 0(1/ In r/) 



Ea;i=l 



P^^''\l - 2xs) ± Pi;''\l - 2x,] 



+ 0(l/lnr/). (5.23) 



This can be seen from the fact that in the limit Inr/ — »• cxd the 'mass' m entering ( p.l9| ) 
becomes large and the kinetic term irrelevant. As seen from ( |5.23| ), the wave functions 
have a two-particle structure corresponding to the relative motion of the conformal spin 
j = 1 quark with momentum fraction Xi (or 0:3) with positive helicity (the first or the 
third quark, in notations of ( p.21| )) and an effective particle with momentum fraction 



^2 + 3:3 (or Xi + X2) and conformal spin j = 2, which is easily recognized as a scalar 
diquark. This is in striking contrast to the structure of the lowest energy state for the 
A = 3/2 baryons ([4.47|) , and in fact the corresponding wave functions are mutually 
orthogonal in the large— iV limit: 

(<undl<')~l/lnr/. (5.24) 

We will further elaborate on the physical interpretation and consequences of this struc- 
ture in Sect. 6. 

Extension of these results to case e = 1 is nontrivial since the higher order corrections 
in qN become dominant, and presents a typical strong coupling problem. From Fig. ^ 
one observes, however, that the quadratic in e behavior of the bound state energy, ( |5.22| ), 
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is replaced by the linear asymptotics starting already from e ~ 0.3 — 0.4. This suggests 
to study the energy spectrum of Ti.{e) in the large-e limit 

e > 1 (5.25) 

and try to find an approximate value of the mass gap at e = 1 by matching the small-e 
and the large-e expansions. This program is carried out in Sect. 5.3 below. 
For completeness, we quote here the results for very small e 

e<l/lnr7, (5.26) 

for which case a simple perturbation theory is again valid and the energy shifts are 
given (up to 0{e^)) by matrix elements of eV over the eigenstates of the nonperturbed 
Hamiltonian. For the states with 0g = (but g 7^ 0) we find, for low lying levels 

£^^\e) - £3/2 = -2g , S'^-^e) - ^3/2 = 0, (5.27) 

while for the states with 0g = ±27r/3 one gets 

S^-'\e) - ^3/2 = -^9, £^-\e) - ^3/2 = -^9, (5.28) 



respectively. Eq. ( |5.14|) is not applicable for calculating the correction to the ground 



state energy with q = since this state is not degenerate. A direct calculation of the 
matrix element (\l/q=o|V^|^l/q=o) based on the exact solutions ( [4.9| ) gives 

S{e)-Eo = -9. (5.29) 

This relation is exact (to C(e)) whereas Eqs. ( ^.27] ) and ( ^.281) are valid up to 0(1/ In^ h) 
corrections. 

5.3 Large— e expansion 

Assuming e ^ 1 is a large parameter, it is natural to invert the logic which we have 
accepted up to now, and consider H3/2 as a perturbation of the spectrum of the Hamil- 
tonian V defined in Eq. ( p.ll ). 

At the first step, therefore, we have to study the spectrum and the eigenf unctions of V 
itself. Although the Hamiltonian V is not integrable and cannot be diagonalized exactly, 
it can be studied in the large N limit using the techniques developed in Refs. ||3^, ^7\ 



for Hamiltonians of similar form. We find that in the leading large— A^ approximation 
the eigenstates in the lower part of the spectrum are given by a linear combination of 
the states diagonalizing the Casimir operators Lf^ and L23. Their relative coefficient is 
fixed by the requirement for the eigenstates of V to have definite parity with respect to 
Vi3 permutations 

1 



1 
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^f)'(a;i,a;2,X3)±^f^'(xi,X2,X3)J + O(r') (5-30) 

^f ^'(xi, X2, X3) ± ^r^'(x3, X2, xi)] + (9(^2) , (5.31) 
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where ^^ (a;,) are functions of the conformal basis defined in ( p.34| ), which we assume 



here to be normahzed as (\1/^ |\l/„ ) = 6kn- The corresponding eigenvalues, £y 



{'$^^l\eV\'$^^l) can be evaluated using (|3l3^ ) and (pC^ as 

1 X V ^'^.v^^'mi- ^^kk 



e-'s'^\k) 



{k + l){k + 2) 



N 

E 



(m + l)(m + 2) ^ (A; + 1)(A; + 2) 



+ C(r/-^). (5.32) 



Taking into account Eqs. (|A.15|) and ( [A.19|) we find 



£'^\k) 



er/-2 [1 ± (-1)^+^] {2k + 3) + C(r/- 



(5.33) 



(A; + l)(A; + 2) 

Here, the integer A; = , 1 , 2, . . . enumerates the levels of V with definite parity starting 
from the lowest one. We observe that for given conformal spin A^ the nonleading 0{ri~'^) 
corrections vanish for each second level leading to £y (k) = ~e/[(k+ l)(k + 2)] provided 
that k + N = even (odd) for levels with positive (negative) parity. It can be shown that 
this result is exact to all orders in l/rj. 

The nonleading (9(?7~^) corrections remove degeneracy of the levels 

£^^\k) - £l-\k) = -2e (2A; + 3)(-l)^+V' + Oir]-') . (5.34) 

Thus, in the lowest part of the spectrum, k ^ N , the eigenvalues of the Hamiltonian V 
belong to the two trajectories, 

£v,\ip{k) - 



{k + l){k + 2)' 
£vAo^n{k) = £v,upik) - 2t{2k + 3)r/-2 + 0{r^) 



(5.35) 



parameterized by a nonnegative integer k. Parity of the eigenstates alternates along each 
trajectory. 

At the second step, we evaluate the matrix elements of 7^3/2 over the eigenstates of 
V defined in Eq. ( |5.31| ). Using ( |A.21| ) and ( |5.31| ) one can write 

N 



1 

2 

1 r 

2 



^Sl^3/2|^K) + (^yi^3/2|<i) = e{k) + 2Y,{-^r^'^kme{m)n^,{->M) 



m=0 



N 



['^V,i\'>i^/2m,i) - (^Sl^3/2|^i7i) 



2e{k)Vtkk + ^ Vtkme{m)Vtmk, (5.37) 



m=0 



where e{m) = 2[^(m + 2) — '?/'(2)]] is the two-particle energy. The sums over m are dom- 
inated by contributions of x = m/N = 0{1). Therefore, replacing e{m) = 21n(A^x) — 
2ip{2) and using properties of the Racah 6j— symbols, ( |A.4] ), ( |A.15| ) and ( |A.19|) , we get 
for the parity-averaged spectrum of 7i(e) 

£W(e) = £[^\k) + (vl/g|H3/2|*g) + 0(l/e), (5.38) 

(<il^3/2|<i) = 41nr^ - 6^(2) + 6^(A; + 2) - 4V^(2fc + 4) +-^ +^^ +(^(r'), 
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while the energy sphtting between the eigenstates with opposite parity equals 

£W(e)-^(-)(e)= 4+^(6) -4-)(e) + 

+ 4 [Inr] + 24j{k + 2) - 3^(2)] 7]-\-lf+''{2k + 3){k + l){k + 2) . (5.39) 

For large r] the correction terms in both cases are dominated by the first term ~ Inr/ 
which can be related to the 'ground state' energy Eq of the Hamiltonian 0.3/2 defined in 
( |4.45| ). It provides an overall shift of all levels and can be absorbed into the definition 
of the nonperturbed Hamiltonian V. Then, comparing the matrix elements ( [5 .381 ) with 



the level spacing in the V^-spectrum in (|5.33| ), we conclude that the large— e expansion 
is well-defined for e ^ 1. 

In particular, for the lowest level A; = we obtain: 

A(6) = ^(6) -E, = -^-e+^ + 0{e-'). (5.40) 

One should stress that all the above expressions are only valid for the lowest levels k <^ N 
in the spectrum. 

It is instructive to examine the flow of the energy levels defined by the perturbative 
expressions ( |5.33| ) and ( ^.38|) . Varying e from e = 00 towards e = 1 we find that the 



energy levels are changing linearly in e, with the slope depending on integer k and on 
the parity of the level. The perturbative correction generates the e— independent shift of 
the trajectories whose amount again depends on the parity and k. The 'critical' values 
of e at which the hnear trajectories cross the ground state energy Eq = £^3/2(0' = 0) set 
up the low boundary for e such that the 1/e— expansion is applicable. Using ( ^.33|) and 
( ^.38| ) we find the corresponding intersection points as 



ecrit(fc = 0) = 0.66, e„it(fc = l) = 5.60, e^^-.tik = 2) = 1Q.97 . (5.41) 

These values are in a good agreement with the numerical solutions shown in Fig. ^. 

Finally, notice that the 1/e— expansion can also be applied to describe the flow of 
energy levels for e < but its range of the applicability is in this case only — 00 ^ e ^ 

— Inr]. The main difference between positive and negative e is that in the latter case 
the lowest energy levels of the Hamiltonian 7i{e) are rapidly approaching each other for 

— Inr/ <^ e < 0. As a consequence, the naive 1/e— expansion becomes divergent due 
to small denominators and should be replaced by the so-called 'degenerate perturbative 
expansion'. 

5.4 Estimate of the mass gap 

Having derived the asymptotic expressions for the mass gap A(e) in the two limits e -C 1 
and e 3> 1, Eqs. ( p. 22] ) and ( |5.4U| ), respectively, we can make an estimate for e ~ 1 by 



matching the two expansions. To this end, we design an interpolating Fade formula: 



e^(fi; + ae) 



^^^ l + 2(/t + 2a/3)e + 2ae2' ^ ' ^ 
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with K = 7r^/(72C(3)) = 1.12549, which reproduces ( |5.22| ) and ( |5.40| ) in the appropriate 
hmits and contains one free parameter a that has to be positive in order to avoid spurious 
singularities. Allowing a to vary within the extreme limits < a < cxo we get 



A 



1/2 



A( 



1^ 



-(0.30000 - 0.34620) 



(5.43) 



which compares very well with the result of the direct numerical calculation of the parity 
average energy of the lowest eigenstates at A^ = 300: 



c{+) 
^1/2 



^1/2 



— Eq 



-0.32097. 



(5.44) 



The results of the numerical calculation of the lowest few eigenvalues of 7ii/2 are 
shown in Fig. [1^. It is seen clearly that the distance between the two lowest eigenvalues 




Figure 10: The difference £^1/2 — Eq for the first few energy levels. 

and the 'vacuum energy' Eq approaches a constant ~ 0.3, while for higher levels this 
distance decreases as 1/ln^iV. Notice, however, that the distance to Eq and the level 
splittings for higher levels are still quite large for InA^ ~ 5 (A^ ~ 10^). The reason for 
this is that the expansion parameter for the upper part of the spectrum proves to be 
1/lnA^ (rather than 1/A^) and the asymptotic large— A^ limit is, therefore, approached 
very slowly. 



6 Distribution amplitudes 
of results 



3/2 

A " 



1/2 
A ' 



1/2. 

A^ • 



Summary 



In this section we give a short summary of the results of phenomenological relevance for 
the physical baryon distribution amphtudes defined in Sect. 2.1 and discuss an overall 
physical picture that emerges. 
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The Brodsky-Lepage evolution equation for the hehcity-3/2 distribution amphtude 

3/2 

(j)^ is exactly integrable and is considered in much detail in Sect. 4. The physical 
interpretation of integrability is that we are able to identify a new 'hidden' quantum 
number which distinguishes components in the A-resonance with different scale depen- 
dence. The anomalous dimensions and the eigenfunctions can be calculated in this case 
exactly using a simple three-term recurrence relation given in Eq. ( [4.3| ). The coeffi- 
cients Un define the expansion coefficients ( |3.32| ) for the eigenfunctions of the evolution 
equation over the complete set of mutually orthogonal conformal polynomials ( |3.42| ) and 
the corresponding anomalous dimensions are given in terms of the same coefficients by 
Eqs. ( [4.43[) and (|2.33|) . Alternatively, we have derived a systematic WKB-type expansion 
for large values of A^ which provides one with a systematic expansion of the eigenvalues, 
[see Eqs. (F^ ), (g), (|i§, { W^ \ and the eigenfunctions [see Eqs. ( ^TTU^ ), (lOUTI ), 



(|A0|)], in powers of 1/A^. 

The case of 0^ is still specific as compared to the general treatment in Sect. 4 in 
that neglecting tiny 5'f/(2)-flavor violation effects due to quark masses (p^ is totally 
symmetric in all three arguments. As a consequence, only one third of the existing 
multiplicatively renormalisable operators have nonvanishing matrix elements, namely, 
those corresponding to the unity eigenvalue ^ = 1 of the cyclic permutation operator. 
Note that the value of 9 alternates along the trajectories shown in Fig. so that each 
third of the eigenvalues gives a relevant contribution. 

Recall that each eigenvalue in Fig. ^ (except for the lowest one for each A^) is double 
degenerate. The two degenerate eigenstates can be chosen either as eigenstates of the 
Q-operator with opposite sign eigenvalues q and — g, or as states with definite parity, 
defined in Eq. ([4.201). The latter choice is more convenient since the parity eigenfunctions 

3/2 

are real and contributions to (j)^ of the operators with negative parity vanish identically. 
One is left with the sum over positive parity eigenstates with real coefficients. 

The most interesting result concerns the structure of the eigenstates with the low- 
est eigenvalue (anomalous dimension) for each A^ which present, therefore, the leading 
contributions to the distribution amplitude in the formal /i^ — *> oo limit: 

XiX2X^^!%'^{Xi,X2,X^) = (6.1) 

= xi{l - xi) Cy,l-^{1 -2xi) +X2(1 - X2) C^Jrli{l - 2x2) + x-i{l - x^) C^Jfl^{l -2x3), 

see Fig. |lT[ The corresponding eigenvalues are known exactly and are given in Eq. ( |4.45| ). 

The physical interpretation of such 'ground states' is most transparent in coordinate 

space. Neglecting the operators with total derivatives, which amounts to going over from 

A]) to the distribution function ^ of one variable in Eq. ( [4.47| ), one can represent the 



three-quark 'ground state' in a concise form as the nonlocal light-cone operator [p8| 

1 /"^ 

^3/r°^(^i' ^2, ^3) = 2 XI ^'^^ I ^^ ^^I(^«^) ^?j(b^« + (1 - v)zp\n) iqlizpn) (6.2) 



ci,/3=l,2,3 
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N=2 



N=4 




Figure 11: Contributions to the A = 3/2 distribution amplitude (j)^ (xj) with lowest anomalous 
dimensions for N = 2 and A^ = 4. The normalization is arbitrary. 



The Tailor expansion of the forward matrix elements of ( |6.2|) at short distances, 2:12, -^32 — >■ 
0, generates the series of local multiplicatively renormalizable three-quark local operators 
with the lowest anomalous dimension for each even N 



-B3/2 {Zi,Z2,Z3) 



A''=even 



'12 



~r ^23 ~r ^^ 



(N + l) 



N 
31 ^ 



Nl=oidu 92,83] 



B{zi,Z2,Z3) 



Za—O , 
81+82+93=0 



(6.3) 



Note integration in (|6.2|) with unit weight over the position of the quark in the middle 
that goes in between the light-cone positions of the other two quarks, up to permutations. 
If renormalization of the operator is interpreted as interaction, integration with the unit 
weight can in turn be interpreted as the statement that the quark in the middle is 
effectively 'free': In the 'ground state' with the lowest 'energy', the interaction of the 
quark in the middle with its right and left neighbours exactly compensate each other. 



1/2 



'1/2 



The evolution equations for helicity-1/2 distribution amplitudes 0^ and 0^ differ 

3/2 

from the evolution equation for (^^ by the additional contribution of gluon exchange 
between the quarks with opposite helicity, see Eq. ( |2.26| ) and Fig. 0. The added terms 
destroy exact integrability, but, as we found, can be considered as a small perturbation 
for the upper part of the spectrum. As a consequence, there is a direct correspon- 
dence between eigenoperators and anomalous dimensions for helicity-3/2 and helicity- 
1/2 distributions and the corrections can be calculated to \/N^ accuracy using the stan- 
dard quantum- mechanical perturbation theory, see Eq. ( p.l2| ). The splitting between 
the eigenstates with opposite parity proves to be exponentially small at large N , see 
Eq. (|TlD. 

For low- lying levels the situation turns out to be dramatically different. We find 
that the two lowest eigenvalues (anomalous dimensions) decouple from the rest of the 
spectrum and in the limit \nN -^ 00 are separated from the other eigenvalues by a finite 
constant A ~ —0.3 ( p.43[ ) that we call the 'mass gap'. The corresponding contributions 
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to the distribution amplitudes are given by 



(6.4) 



,(3,1) 



where Pjy ' are Jacobi polynomials, and correspond, in the same sense as above, to the 
contribution of the nonlocal light-cone operator 

5(^1,^2,2:3) = 



Jjk 



iyiqj flq]){zin) fiql{z3n)6{z2 - Zi)± fiqj {ziu) {yiq\) fiql){z3n)6{z2 - Z3) (6.5) 



Formation of the mass gap in the spectrum of anomalous dimensions is, therefore, nat- 
urally interpreted as due to binding of the quarks with opposite helicity and forming 
scalar diquarks. 

Note that while the expression for the eigenfunction in (|6.1j ) is exact, the result in 
( |6.4|) is only valid in the asymptotic InA^ — i> 00 limit. In the coordinate space picture, 
the restriction to large A^ is translated to the condition that the light-cone separation 
between the same helicity quarks is very large to allow for the formation of a diquark. In 
momentum space, the result means that at sufficiently large scales Q^ the quark carrying 
a very large momentum fraction is more often with the same helicity as of the parent 
baryon. This observation seems to be in qualitative agreement with phenomenological 
models of baryon distribution amplitudes derived from QCD sum rules [|p, |2( 






Figure 12: The 'ground state' eigenfunctions ^(x) for the Hamiltonians TC1/2 (solid), 7^3/2 
(long dashes) and V (short dashes) for N = 19. The normalization is to unit integral 
Jdx^{x) = 1 and to the unit first moment fdxx^{x) = 1 for the symmetric and the an- 
tisymmetric wave functions, respectively. 

One has to keep in mind, however, that the diquark picture of the states with the 
lowest anomalous dimensions only becomes quantitative for very large A^. To illustrate 
this point, we plot in Fig. |1^ the exact eigenfunctions \l/ for 7ii/2 at A^ = 19 (with 
positive and negative parity) corresponding to the lowest anomalous dimensions (solid 
curves), and compare them both with the lowest-level eigenfunctions for 0.3/2 (dashes) 
and with the 'diquark' eigenfunction corresponding to Eq. (|6^ ) (dots). It is seen that 
the eigenfunctions of ?^i/2 for this value of N are still very close to 7Y3/2 and only start 
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to develop small 'horns' close to the end points, characteristic for the diquark picture. 
As mentioned above, the true large- A^ limit is approached very slowly since parameter 
of the expansion is in this case 1/lnA^ rather than 1/A^. 

7 Conclusions 

To summarize, in this paper we have developed a new theoretical framework for the 
description of baryon distribution amplitudes in QCD, based on integrability of the 
helicity-3/2 Brodsky-Lepage evolution equation. The mathematical structure of the 
evolution equations reflects a clear physical structure of the distribution amplitudes that 
we tried to emphasize. A lot of analytic results is obtained, in different limits. 

The formalism proposed in this paper is rather general and can be applied, as indi- 



cated in |T^, to the studies of quark-antiquark-gluon and, possibly, three-gluon distri- 
butions. 

Three general questions related to the theory of three-particle distributions are not 
covered in this work and deserve further attention. First, as we have indicated, ana- 
lytic continuation of the spectrum of anomalous dimensions of three-particle operators 
to the complex angular momentum plane is intrinsically ambiguous. One has to study 
whether this mathematical ambiguity is resolved by imposing certain physical conditions 
on the amplitudes. Second, the solution of the evolution equations for three-particle dis- 
tributions depends on the nonpertrubative initial conditions. Depending on the choice 
of the three-particle distribution amplitudes at low scales, there is a possibihty that at 
large evolution times the solutions to the evolution equation become independent on the 
initial conditions and are governed entirely by perturbative evolution. Such perturba- 
tively driven distribution amplitudes would generalize the GRV partonic distributions 
Il39| which prove to be successful in the phenomenology of hard inclusive processes. Fi- 
nally, the integrability of the evolution equations reveals an additional hidden symmetry 
of QCD and its close relation to exactly solvable statistical models. Remarkably enough 
the same symmetry has been observed in the studies of the Regge asymptotics of three- 
particles distributions. These properties are not seen at the level of the QCD Lagrangian 
and their origin needs to be understood better. 
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A Appendix: Racah 6j-symbols of the SL{2) group 

The functions \I/)vn i^i) introduced in Sect. 3.4 have a simple group theory interpreta- 
tion. They define the addition rules for the sum of three conformal spins j = 1 each 
corresponding to the SL{2) generators La^k (^ = 1,2,3) acting on the light-cone coor- 
dinates of three quarks and the superscript (12)3 indicates the order in which the tensor 
product of three SL{2) representations has been decomposed into the irreducible compo- 
nents for which '^j^^ (^«) is the highest weight. According to ( p.34| ) and ( |3.36|) , ^j^n (^0 
describes the irreducible component for which the total conformal spin is h = N + 3 and 
the conformal spin in the channel (12) is equal to ji2 = n + 2. Changing the order in 
which the spins are added one obtains an equivalent basis of functions \&„ i^i) that 
are linear related to \l/„ (xj) through the Racah decomposition ( |3.37| ). The coefficients 
Qkn thus define the 6j— symbols for the discrete positive series of the SL{2) group. It 
is easy to see that the corresponding basis functions are transformed into each other by 
cyclic permutations 

^f!^\x^,X2,X^) = V&iif (X2,X3,X1) (A.l) 

and, therefore, the (A^ + 1) x (A^ + 1) matrix Vt represents the operator of cyclic permu- 
tations V ([4.1 0|) in the conformal basis 



(V^(12)3|^l(23)^ = {^'^'''\V\^T'') = ^nmj^^y (A.2) 

where we used the expression for the norm of the basis vectors (|3.39|) . Since \l/n is a 
real function of Xi, the calculation of the scalar product leads to real matrix elements 

N 
n=0 

where the second relation is the unitarity condition. Since V^ = 1 and V^ = Pi2^^i2 = 
V^^, the matrix Q has to satisfy the following conditions: 

O^ — A O^ — r_Ti"+"*0 — o~^ ( \ A\ 

where in the last relation we used the identity 

Vu'^'jif (xi, X2, xs) = ^]if (X2, xi, xs) = i-lT'^'jif {x^,X2, xs) . (A.5) 

Explicit expressions for the matrix elements Qkn can be obtained in terms of the 
4F3— hypergeometric series of unit argument |40|. To show this, consider the defining 
relation ( p.37| ) and choose Xi = —X2 = x and X3 = 1 so that xi + X2 + X3 = 1 and 
Xi + X2 = and \1'„ (xj) reduces to x" up to a numerical factor. Rewriting the 
Gegenbauer and Jacobi polynomials in terms of hypergeometric functions, one brings 
Eq. ( |3.37| ) to the form 

^/— m, — m — 1 \ ^ /— A^ + m, A^ + m + 5 \ v^ ^ r, /. x 

2Fii '2 xUFii ^^, + 4 xjp„ = 5^fi^„xV„ (A.6) 



n=0 
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where 

CTn = (-1) 



^_„,(iV + n + 4)! 



{n\Y{N - n)\ {2n + 2){2n + 3) ' 



1 ( Kf A_ jj, \ A)\ 



The two hypergeometric functions entering ( |A.6|) are polynomials in x of degree m and 



N — m, respectively. Their product defines the generating function for the matrix ele- 
ments ^mn- Expanding the l.h.s. of [A.6| in powers of x, one can write the coefficient of 



x^ as the 4F3— function of unit argument and identify it with fimn- Explicit expressions 
for arbitrary A^ and n are rather cumbersome, see, e.g., ||41|| . 

We are able, however, to find a simple approximate expressions for Qmn which are 
valid in the WKB limit of large spins A^. To this end, notice that the matrix elements 
Qnm satisfy second-order finite difference equations. One finds them by applying the 
operator L23 to the both sides of (|3.37|) and taking into account ( p.40|) 



n + 2 [n + lj[n + 2j n + 1 



where m,n = 0, . . . ,N and Q^ -i = 0. Let 



n^n = (-1)"^ (A.9) 

In 



SO that Eqs. ( [A.8| ) and ( [A. 3D are replaced by 

2(m + 2)(m + l)ujm{n) = -^{uJm{n + 1) - u;,„(^)) + ;^(u;^(n) - u;„(ra - 1)), 

Y.n=0 j;^m{n)uJl{n) = fmSml, (A. 10) 

ujmin) = uJnim), 

defining the system of orthogonal polynomials ujm{n) in the discrete variable n with 
m = 0, . . . ,N. The initial condition for the recurrence relations in ([A.10|) can be found 



from the relation $^ (zi) = ^^=0^'="'^" (-^*) ^^ ^^^ limit zi — Z2 ^^ using ( p.55| ) 

and (|3.56| ) as 

2('-l)^ 
a;„(0) = ^-—^ -fJ2n + 3) . (A.ll) 

There are two limiting cases when the recursion relations for uJm{n) can be solved 
analytically. In the first case, in the limit 

Tl 

X = — = fixed , TJi = fixed , N ^>- 00 (A. 12) 



the recurrence relation in ( [A.10| ) is reduced to the differential equation on the function 

/, / N 1 / / N 4(m + 2)(m+l) , , ,. x 

X 1 — x^ 
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Picking up the polynomial solution we obtain 

uj^(x) = N^x\l - x^)C^\l - 2x^) . (A.14) 



The normalization is fixed by the second relation in (|A. 10|) where one replaces ^ 



N Jg dx. Finally, substituting (|A.14|) into ( |A.9D we get 



^mn = ^^xC^/2(l - 2x^) , n = xN, (A.15) 



which is valid to 0{1/N) accuracy in the limit specified in ( |A.12| ). 
In the second case, in the limit 

n,m = fixed, N ^ oo (A. 16) 

the recurrence relation in ( [A .101) is reduced to the condition 



(ujmin + 1) - uj„Xn)) = — —-{uJrrXn) - ujrnin - 1)) . (A. 17) 



n + 2 n + 1 

Its solution satysfying ([A.ll|) is given by 



uJra{n) = ^{-lf{n + l){n + 2)(m + l)(m + 2) x [l + 0{1/N^)] (A.18) 

leading to 

Q^^ = A^-2(-i)^+"(2n + 3)(m + l)(m + 2) x [l + 0{1/N^)] . (A.19) 

Having defined the matrix flnm it becomes straightforward to calculate matrix ele- 
ments of the Hamiltonians 7^3/2 and 7ii/2 in the conformal basis. To this end we write 

n-i/2 = Hu + VHuV"' + V^Hi2V-^ (A.20) 

and similar for 7ii/2- Then, applying this identity to a basis function Wn {xi) and 
using the properties ( |A.4| ) of the fi— matrix we find that 7^3/2 can be represented in the 
conformal basis by the {N + 1) x [N + 1) matrix [7^3/2] nfc 

^3/2^^)=^ = Jyiz/2U^T^\ (A.21) 



fc=0 

N 

'I 1 \mo _/ \c 

''mki 



\n^l2U = e{n)6nk + [(-!)"■ + (-1)'] 5^(-l)'"l^„„£(m)fi, 



m=0 



where e{n) is the energy of the two-particle Hamiltonian defined in (|4.40| ). It is easy 



to see that this derivation relies on the two-particle structure of the Hamiltonian only 
and is not sensitive to integrability properties. In particular, the similar representation 
holds for the Hamiltonian 7^i/2, with e{m) in the sum over m shifted by the exchange 
interaction term l/[(r7i + 2){m + 1)]. 
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B Appendix: The low-energy effective Hamiltonian 

for n{e) 



The eigenproblem for the matrix (|5.19|) takes a well known form once we interpret the 
integer /c as a discrete momentum variable. Denoting the corresponding eigenvector as 
c = {cfc} we construct the wave function in the configuration x— space as 

h——k 

n, — "rnax 



The restriction to \k\ < A;max serves to remind that (|5.19| ) presents an approximation to 



the Hamiltonian 7i(e) which is only valid for \k\ -C Inrj. It is natural to expect that 
the lowest energy levels of the matrix ( |5.19| ) are not sensitive to the UV cut-off fcmax, 
or, equivalently, the corresponding eigenstates x{^) ^-^e smooth functions of x at short 
distances Ax ~ l/lni]. To the extent that this 'decoupling' property holds true, which 
we are going to verify a posteriori, the low-lying levels of the Hamiltonian 7i(e) coincide 
with the lowest eigenstates of (|5.19|) so that the latter can be considered as the effective 



low-energy Hamiltonian for the former. Having this in mind, we temporally send the 
UV cutoff A;max to infinity and assume the matrix (|5.19|) to be of infinite size. 



Using the transformation ( [B.l| ) one can map the eigenproblem for the matrix ( p.l9| ) 
into a one-dimensional Schrodinger equation for the wave function x{^)- It follows from 
( P.l|) that x(a;) is a (anti)periodic function of x with the period 3 for even (odd) values 
of N, respectively: 

X{x + l) = {-lfx{x). (B.2) 

The two cases should, therefore, be treated separately. Let us first consider the case of 
even N and split the wave function into the sum of three terms 

Xix) = x+ix)+X-ix) + Xoix), (B.3) 



Xa / ^ Cs/c-j-Q, 



2i7r, 



^^-x{3k+a) 



fc=— oo 



where a = ±1 , 0. Each component presents a (quasiperiodic) Bloch-Floquet wave func- 
tion with the period 1 and the quasimomentum 27ia/3: 

Xo{x + 1) = Xo{x) , x±{x + 1) = e^^'^/^x±{x). (B.4) 

It is straightforward to show that the eigenvalue problem for matrix ( |5.19| ) is equivalent to 
the Schrodinger equation for the three Bloch-Floquet waves Xa{x) propagating through 
a periodic array of (5— function potentials and interacting with each other: 



1 1 °° 

fc=— oo 



^Xa{x) - Yl X(3{x)e 

/3=±1,0 
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^xia-P) 



Exaix) (B.5) 



with 



S{e)-Eo = gE. 



(B.6) 



This Schrodinger equation generahzes the famous Kronig-Penney model of a single par- 
ticle in a periodic 5— function potential and its solution follows the same procedure [^ . 
Namely, the solution to (|B.5|) on the intervals of periodicity n < x < n + 1 with n integer 



are given by the plane waves 

Xa{x) 



ar,[n) e 



2ipx 



+ bJn)e 



-2ipx 



(B.7) 



with the coefficients aa{n) and ba{n) depending on n. The corresponding values of the 
energy take the simple form 

2c/p2 



£{e)=Eo + gE = Eo + 



Eo + 



m 



2C(3) 
In (r^e 



i7E' 



:P 



(B.8) 



The possible values of the momentum p are restricted by the quantization conditions 
that one establishes by requiring Xa{x) to be a continuous function of x satisfying the 
periodicity condition ( |B.4|) and its derivative d^Xai^) to have a discontinuity at x = n 
which can be calculated by integrating the both sides of 



dxXaix] 



n+S 



n—d 



-m 



3xa{n) 



/3=±1,0 



n(a-/3) 



(B.9) 



with 6 —>■ 0. We find that for even N the quantized q have to satisfy one of the following 
three conditions 



smp 

2 , 



p cotp (3 — tan p) 
3 



0, 
3m. 



p tanp 



1 



tan^p 
tan^p 



3 

--m , 
2 



(B.IO) 
(B.ll) 

(B.12) 



which define three different branches for the dependence p = p{m). For odd N the same 
conditions look like 



cosp = , 
p tanp (3 — cot^ p) = —3m , 



p cotp 



3 — cot^ p 
1 — cot^ p 



-m . 



(B.13) 
(B.14) 

(B.15) 



The solutions to (p.lO| ) and ( p.l3| ) do not depend on the perturbation e and the cor- 
responding energy levels coincide with the levels of the unperturbed Hamiltonian. This 
happens because for these values of the momentum p the wave function vanishes at in- 
teger points Xa{n) = and, as a consequence, the interaction term in (|B.5|) vanishes as 
well. 



69 






=, ^_==^^ 


6.0 







=i 


_ 


^ 


_ 




^ 

















"- 


r 


I 


-_ 


_ 


— 


I 


— 


— 


— 


=— 










T~~ 











3.0 


s= 


— 


_ 


._ 




___ 




_ 


_ 


■= 


=— 




"~ 


_ 






— 


— 


-. 


.- 


— 


- 


— 


- 


kc"~ 







0.0 




— 


— 






— 




— 


— 




^^-a 


■K 




























































\ 




3.0 


























s 


\. 



-8.0 



-4.0 



0.0 



4.0 



8.0 



Figure 13: The flow of the quantized momenta p = p{m) for even (left) and odd (right) 
values of N . Dashed, long-dashed and solid lines are described by three branches defined in 
Eqs. ( plCp , ( pup , ( pl2|) and Eqs. (^IsD , ( pip , ([bIsD , respectively. Negative values of \p\ 
correspond to pure imaginary momenta p. 



Since the matrix ( ^.19|) is herinitian, its eigenvalues ought to be reaL Then, it fol- 
lows from ( p.8|) that quantized p can take either real or pure imaginary values. In the 
latter case, the energy E becomes negative and the wave function ( |B.7D is exponentially 
decreasing with x, indicating formation of a bound state. We will see that these bound 
states are precisely the ones that are responsible for the decoupling of the pair of lowest 
levels in Fig. |[ 

To this end, we solve the quantization conditions (p. 11 ) and ( p. 14 ) for different values 
of the 'mass' m. The resulting dependence p = p{m) defined by the three quantization 
conditions in ( [B.ll ) and ( B.14 ) is shown in Fig. |13|. The following comments are in 
order. 

At m = the solutions to ( |B.11| ) and ( [B.14| ) are given by pk = vr(A;/3 + 5n 1'^) with 
k integer and the corresponding energies (p.8| ) coincide with the energy levels ([4.91| ) of 
the Hamiltonian 'Hz/2- All solutions except the one with fc = and A^ = even are double 
degenerate. 

For small m <^ 1, or equivalently elnr/ <C 1, the degeneracy is removed and each 
branch, pk = p\. (m) and p = p\ {m), evolves independently. The slope of the trajec- 
tories at m = can be found from Eqs.( |B.l"OD - ( |B.15| ) as 



d 



P3k±l 

dm 



3 



P3k±l 

dm 



1 

1' 
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d 



dm 



d 



-1 



dm 



n2 



0, 



with p): 



(±), 



m, 



0) = n{k/3 + 57v/2) ^ 0. For even N and m 



non-degenerate level po = which evolves as 

dpi ^ 
dm, 



1 
2 



(B.16) 
we have a single 

(B.17) 



It is easy to see that Eqs. (|B.16|) , (|B.17|) and (p.8| ) are equivalent to (|5.27|) , (|5.28| ) and 
( ^■29| ). Since Po(0) = 0, it follows from ([B.17|) that po becomes pure imaginary for an 
arbitrary small m > and the corresponding eigenstate describes a bound state with 
the energy 

4ound(e) = i^o - ^1^^^^^. (B.18) 

with e ^ l/lnrj. For even A^ there exists the second bound state which is formed for 
nonzero mass m, > rricrit- The value merit corresponds to the nontrivial solution of ( |B.11| ) 
at p = 



TT 



merit 



9C(3) 



ln(?7e'^^) = 1 . 



(B.19) 



The similar phenomenon occurs for odd N. In this case, two bound states are formed 
for m > and the corresponding critical values of the masses rricrit (or equivalently ecrit) 
can be found from ( [B.14 ) and ( B.15 ) for p = a.s m,crit = 1/3 and 2/3, respectively. It 
is easy to see that for m close to merit the mass gap, inbound (e) — Eq depends linearly on 
the perturbation e for the both bound states. 

To understand what happens with the spectrum of the Hamiltonian as e varies, 
consider the solutions to Eqs. ( p.lO|) - (p.l5| ) in the two extreme limits: m = — oo and 
oo. The real quantized pk are the same in both limits. 



Pk 



-nk 
4 



k = 1,2 



(B.20) 



and, as a consequence, the energy levels in the 'continuum' gE > are also the same and 
are given by (|5.21| ). Moreover, for large positive m there are additional pure imaginary 
solutions to Eqs. ( |BTTD , ([KT^ ) and Eqs. (|KT^ ), ( |B7[5|) 



TT 



^Pbound 



4 



-m 



12C(3) 



ln(?7e 



,7-B^ 



(B.21) 



which give rise to two bound states with the energy given by ( p.22| ). 

Thus, the flow of the energy levels from m = — cxd to m = cxd is such that the 
continuum stays unchanged and two lowest levels of the continuum 'dive' into the vac- 
uum. These two bound states become separated from the continuum by the mass gap, 
Abound (e) — Eq, whosc sizc grows linearly with e at small m ~ merit and quadratically at 
large m. 
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Recall, now, that Eq. ( |5.19| ) presents a low-energy approximation to the Hamiltonian 



Ti.{e) and we have to check that values of qN are small on the solutions. Since, according 
to (|5.18|) , the kinetic energy contribution to the Schrodinger equation is O {{qN)'^), the 
condition qN < 1 can be expressed as 



/ dx 

6m Jq 



9M^) =£E^'N'= ^ <!' (B-22) 



2771 



2g^ 
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cf. Eq. ( [B.8|) , with the eigenstate xi.^) normalized as |/q dx\x{x)\'^ = 1. Since 771/ g 
scales at large A^ as ~ In 77, the restriction in ( |B.22| ) imposes the UV cut-off on the 
quantized momenta of the states \p\ < In?]. It follows from ( B.20|) that for the states 



in the continuum, gE > 0, this condition is satisfied for the k = 0{ln7]) lowest states 
only. For the two bound states with gE < the relations ( P.22|) and ( p.21|) lead to the 
condition 

^m ~ e^ < 1 . (B.23) 

Thus, at large N, our assumption about decoupling of the low-energy levels (smallness of 
qN) is justified provided that e < 1. For e ~ 1 the higher-order qN corrections to both 
the kinetic energy and the potential terms in (p.5| ) become significant and eventually 
start to play the dominant role for e ^ 1. In this case, as it was shown in Sect. 5.3, 
the eigenstates of the exchange interaction V provide the appropriate basis and the 
Hamiltonian 7^3/2 can be treated as a perturbation. 
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